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TI. INTROTUCTICON 


mr 6CONCEPT 

The intent of this thesis was to agveiop a computer 
meoeram that en experienced controls engineer could use _ to 
Metro him mod€l, simulate, and realize contrcl systems with 
Mmemeconvenience of 2 aesatoz microrcmouter. Tre cortrol 
engineer Vil eNGt Needs to be an experiansei computer 
m@moprammer. In order to accomplish the above requirements it 


is necessary tnat the trogram be easy to use, effectively 


by 


Mecuiring that tae program write tae computer cede itsel 


Mmeeelittie help from the control designer. This is cemmonly 


ws 
+ 


Mererred to aS automatic programmines. i121 order 


tU 


MeOSTaM to v2 4 useable tooi for the controls d2signer the 


cr 
t- 
Ut 


motem must also produce a good motel of tre system tra 


5 ad 
0) 
‘D 


Meine designed. Cbviously all control systems cannot 
meesidsred in this tyve of project. Tre anproack téaxer iin 
mrs thesis "was to use control systens trat rave 2rror 
Metals designed to be eaual to zero or differ fror zero by 4 


Mumeemoanit. OLner sidelines of taHis project include: tine 
Peeevions must be simple difference enuations, the syster 
meee ailow for calculation delays within tre computer and 
Meese delays must be corrected ror automaticaily or at isesit 
mmeewe fOr the possisility that the designer can correct for 


meemaelays, and the rrogram must te able to medel tre plent 


te) 





between the sample perioas of the computer which tasically 
means there are two sample periods, one for tne somputer and 
one for the plant. Since the plant cannot be programmec as 
mal og it must be discretized. Using a faster sample rate for 
the plant allows simpler equations yet mainteins a good 
Mmemresentation for the plant. additionally, the program must 
be able to model analog to digital and dizital to analog 
conversions. After the model ras teen designed and the 
Semilation run, the results that are available from the 
program must be easy to use anc eesy to interpret by the 
mootrols engineer. In eard]r to meet ail Oe tnese 
requirements, most of the program devslopment itseif ere pe 
Mransbarent to the controls engineer. Tie arproack ta 
meme the control engineer to orovide oniy equations. funciton 
mames, and function relationships by providing tne names * 


@eger Iunctions used aS invuts 


Pee eooUMPTIONS AND RECUIRZEMENTS 
mie assumptions and requirements are siated below. 
aS Mentioned above, tne program must se reiatively easy to 


Meer rce voroeram must have high precision math capatilit; 


ht 


pmece it was aesiened to be used with tontrol systers. 37eec 
Mmemecompliling and oneration was not considered critical at 
mmo point since Micropro:essors are hbesoming small2r aac 
faster which will solve any reelization problems in the aear 
Mmeurs. A Major consideration is that ence *ne sintulations 


Bee been run and the results are acceptatl? to the designer 





then the simulation equations should be easily converted into 
realization equations. Basically this means the simulation 
equations in the computer should be easily converted into 
hardware equations that can be used to control areal plant. 
Me order to accomplish this taSk it was decided that the 
program must develop source code that could be compiled ty a 
high level language compiler for either computer simulation 


or hardware implementation. This approach makes the prerram 


Meansportable and hardware independent. 


ee 6AM PLEMENTATION 


irom ore ject waS Implemented entire GN Ml 210 pee cess OF 


os 


Mmesed systems to show that it 1s feasible ana practical to 


O 


use these small systems for major development. The tresent 
mesavem iS slightly constrained due to the leneth of sompile 
time. All programs were developed usine a Z&¢ microprocesscer; 
mowever, the code is written in high level source lanesuvaces 
therefore the programs will also run on the simpler E¢&¢ 
Mmaseac microcomputers or on the more complex sixteen bdii 
memeroprocessors such as the 89EE. 

The ianevages chosen were PL/1-&@ from Ligitel Research 
Because of the math and hardware cavatilities «f the language 
mmeewines for control of analog to digital and other processes 
Merde the computer, and DBASzt II from Ashton Tate whith 
a@ilows easy and high level manipulation of databases and file 
memmems., One of the major points of this tnesis 15 to 


Manipulate fume t ons and functional relationshins as 


a 





databases. Once the relationships have been developed, the 
database program writes PL/1 source code from the database 
information. The PL/1i source code is then compiled ty the 
PL/1 Compiler and the simulation run on the resident 
hardware. 

Pveeuresent Drograms and the test example in this report 
provide simulation only and send the results in floating 
mene, format to a disk file. The floating point file can then 
Mmemeconverted to display the information in tabular or 
Zraphical form such as shown at the bDacz of this report. The 
information can be displayed in a variety of forms suck as on 
msm lay consoles, printers, plotters, or graphic terminals. 
moe real power for design is in the use of graphical 
equipment to provide the control engineer with tne entire 
picture of the system under test. The Hewleét Packard ¢8&72B 
Meeotter and 2<647A sraphics terminal and the Intercolor 8363 
meeipnics terminal were used to snow that the data produced 
was basically machine independent. The code can he compiled, 


Mem, and displayed on a variety of different machines. 


t— 
A] 





II. GENERAL DISCUSSION 


A. CONTROL ENGINEER INFORMATION 

The control system engineer does not need to know how the 
Bemputer programs work in order to use the system. it is 
assumed that the control engineer is a novice prosrammer with 
some Knowledge of the CP/M operating system and tlimited 
Beckeround in program writing; however, it is essumed he is a 
Knowledgeable control designer. This means te is fairly well 
versed in taking control systems and producing discrete 
System equations from the continuous control equations. These 
Squdations are designed and optimized by whatever means are 
available either by hand or on a mainframe computer. 
meaorocomputers Will soon be able to ‘to major system 
development equations.) Once the equations are developed the 
Mesiener enters them into a microcomputer and tests the 
mmeations. [f the simulations are acceptacle then it srouid 
Memeeeah §86Fasy matter to convert the simulation code for tke 
Sontroller ees en into development code meeve ant ia 
implementation into rardware. 

For ease of use tne drograms are menu driven and use 
maple Cursor arrows or commanas for pointing at selections. 
Additionally very basic information is needed sucn as thea 
mame of functions and e@quations ina simrle format the 


memmputer can use, This procedure is accomplishea by usine tne 


> 
CN 





database Drosram which only allows a user to enter 
moecormation in very Spo Ciebeownieids Of every specific type. 
The type of information required from the designer is to add 
Mee delete equations, to change initial conditions, and _ to 
provide function names. 

There is an initial learning curve with this program but 
Sere attempt has beén made to keep 1t as short as ovossiole. 
There are some peculiarities put in as a requirement to mae 
interfacing between the human anc the computer possibdle. some 
Bmemples include names such as INPUT:A or INPUT:E that 
indicate which function the designer is talking accut to the 
Computer. As an example INPUT:A will associate a functior in 
mee computer with the function name the designer nas selected 
men as INPUT:A = Un. Another very abstract nam> is AiPret. 
Tre user is expected to enter either iif he wants tre 
MesoeNtT value of the function named under INPUT:A or 2 i? hte 
wants the DELAYED value of the function named under [NPUT:A. 
If the desizner will be using both the present and fast 
values then enter 1 to get the present value and the vast 
maeue Will also bé available. The strange requiremert fer 
AlP2D is vrovided only to allow much faster prcecgrams when the 
code is implemented in hardware and the designer only needs 
past or delayed values of a function. Two other codes ar2 IAP 
See tAL for INPUT A PRESENT and INPUT A LTELAYEL. These two 
codes are provided to Keep the equations neeiea for the 


functions short. Another variavle is 0 which is usec to 


14 





define the delayed or previous time sample of the functions 
own OUTPUT. There are similar definitions for up to five 
Momouts INPUT:A through INPUT:E for each function. A function 
can be defined in three ways; °“F means the function is an 
GCquation that the computer finds a register to hold the 
result in, ‘ALT(x) means the value is to be placed in one of 
eight analog to digital registers, and TA(x) means the 
value is to be placed in one of eight digital to analog 
registers. An example equation and simple tlock tiagrar is 
shown in Figure 1. 

Other variables used include “tT” for the discrete sample 
period and T1 for the sample period of the plant. If the 
maeue of T is normalized and Ti has the same value as T then 
mee programs will run faster in the realization of the code 
since these values effectively disappear from the equations. 
It should be noted that all functions are done in floating 
point so the simpler the equations the faster the ccde. The 
example problem at the end of this paver normalizes T but 
meres includes it in the equations to show how it could show 
Meeeein the function. Finally the designer is allowed to 
Mmeeriy End criteria ana intervals for calculating tre plant 
Mmewative to the computer itself. 

After the designer has entered the equations into the 


database and is satisfied that all function relationships are 


Mmmmect then ali that is required is to enter the correct 


a 


Menu selection that compiles the code. The computer progra 





IAP TAD 





Lhe 


NOTES: 

Assume by previous definition that INPUT:A = X, INPUT:5B = ff, 
FUNCNAME = Dz, and the below equation is desired: 

Meme X(k) — (2 ¥°xX(k-1)) + ¥(k) - (@.5 * D(k-1)) 

Making substitutions: 

Meee = Dik), [AF = X = X(k), IAD = X{k-1), IBP = Y = Y{k), 
mage O = OUTPUT TCELAYED = D(k-1) leads to the equation --- 
Meera) - (2 * LAD) + IBP - (8.5 * 0) 


Notice that the above function only haS two inputs, X ens Y 


Figure 1 


Simple Difference Equation [emonstrating System Variables 
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will now write the necessary source code from the database 
and subsequently invoke the PL/1i compiler which will turn the 
code into a program that will run on any of tne previously 
discussed CP/M based systems. The executable file produced is 
called ANSWER.COM and is run on a CP/M based system by typing 


DEO mpc. 


"ANSWER in response to the operating system input 
Puring program ANSWER.COM execution the results will be sent 
mea t:lLoating point rile called DATA.¥LT. Also, while the 
program is running the intermediate results are printed on 
the display device for program troubleshooting should the 
results not te as expected. The file TATA.FLT is not 
printable; however, an example utility program called 
PRINT.PLI is provided that changes the floating point file 
into a printable text file. This program is executed by 
typing PRINT FILE where FILE is the name of the designer 
provided text file the designer wants the results placed 
Into. With further programming development the floating point 
data can be converted into any format needed by different 
Meaeroes SUCH as plotters aS Shown by the plot results in 
Appendix C. 

The information for understanding the database program 
Mmeeration is provided in Chapter V¥. This section is not. 
necessary for understanding at the design level but it does 
eepidain variables that will be needed in equation 
@evelopment. Additionally the example problem will be heipful 


1 
hak Nee 


ct 


to the design engineer. when the designer is ready to use 





program all that is needed is to enter DBASE MENU’ to the 
CP/M operating system prompt. AS previously stated the 
designer will need some familiarity with CP/M before the 


program can be of use. 


B. PROGRAMMER INFORMATION 

This section covers the experience necessary for aaeereln 
maintenance and further development. These programs would be 
worked on and coded by an individual that kas a general 
control system background; however, he must have substantial 
experience with computer programming. The type of computer 
experience needed pertains to operating systems, event 
counters, and automatic programming. The programmer must have 
a good knowledge of both PL/1-8@ and [TEASE II. Additionally 
it would be necessary tO know microccmputer hardware 
implementations such as analog to digital converters and 
Sample ana holds. These devices are simulated in the computer 
programs but further useful extensions would involve 
implementing the code into realized hardware. The information 
required for the programmer to understand the programs is 
Semrained in the remainder of this tnesis. The majcrity of 
meas thesis is the computer code required to implement a 
design environment for a novice programmer through the use of 


dutomatic programming. 


ro 





ipo LOS WNC RAS IES 


—— Pee ae ee oe oa ee Se ees ew ewe See 


Major problems encountered in this thesis had to do with 
meoerdmming and efficiency. The major problem is that there 
is no one programming language that contains all the features 
Mecessary to implement the entire project. JLUBAsS# If was 
chesen because of the relational database characteristics 
even though Coes every eG ed math a6 hardware 
Capabilities. The requirements for reasonable matn is the 
reasor PL/1-58 was chosen and though some database systems 
are avaliable for PL/1 they do not have all the features of 
Merril, A major deficiency of [BASE II is that it will not 
allow GMemEprinting =» Of ia semicolon. ine “sericolon 1s 
considered a command character in [EBASn [I and therefore 
Meservyed. This can be considered a major bug for a database 
eyotem, Papecial yeeome thar was designed for business 
mpoeeycations. It is roolish to s2e€lieve business will not want 
to mc Lude a semicolon someplace in one of Ler 
Mmeamsactions. Ashton Tate was notified and they replied it 
Mero WOTLH trying tO fix. Tnis presents ea problem since 
the databdvase program was used to write the source code in 
Fo/1 and PL/1 reauires a semicolon at the end of every 
Statement. The vroblem was worked around by kaving I[BAS=® ft! 
Write a backslash every place a semicclon was needed in tae 


Mey source text. Next ¢ PL/1 preprocessor called REFCRY.PLI 
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was used to change all backslashes to semicolons prior to 
invoking the PL/1i compiler. This takes a large amount of the 
eeppile time “which could have been easily corrected by a 
reasonable response from Ashton Tate. PL/1 also has some 
problems; however, Tigital kesearch was more than happy to 
meee tO overcome deficiencies. ‘There are several problems in 
their floating point routines due mainly to conversion rules 
meach can cause strange effects in the results. Most floating 
Mernt functions such as the exponential function tend tc 
convert the arguments to fixed before calculation and tkis in 
effect truncates the arguments. AS an examvle if one desires 
meumeraise an exponent to a power and that power is a funstion 
Such as X/Y then both X and Y and the result get truncatiei. 
Mme tks 6alCUcvtraction less than 1.¢ the effect is to truncate 
memes result to €.6 and tne resvit of the exponent to a 
M@eetcioOna! power always endS up “qual to 1.€@. The way around 
mers Droviem is to convert X and Y to CHARACTERS before doing 
Meee aGdivision. This 1s only vart «cf the problem. Even if X and 
meorisionally were ‘floating point and were convertea io 
Memeracte’rs the result of the exponent evaluation is “fixed 
decimal. Now this result must ve converted by anotner 
Mmemevion called #7C in order to get the resuit back to tze 
meovem in floating point. An example of this shows up in the 
Mem TOR.PLI program for the function callea CF. This does 
moe some limitation on the type of efuations that the 


Meeeener Can use without knowing some peculiarities of tke 
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designer can use without knowing some peculiarities of the 
system. 

For the sake of program efficiency, as far as development 
goes, all variables or functions the designer uses are 
considered floating point. This makes programming easier with 
Myener precision; however, the code takes longer to execute 
even though the execution time was not considered critical in 
this implementation. Using floating point for all functions 
introduces some inaccurate realizations when dealing with 
moaloe to digital and digital to analog converters. This is 
not a major problem since these devices will have greater 
meemracy in the future or the devices can b2 moceled more 
@ecurately in future designs. 

meee results presentiy produced by this system are placed 
on floppy disks inside the computer. writing to the disk and 
displaying the results on a terminal are an obvious slowdown 
an System operation. These inefficiencies do not show up in a 
Mearization of the control system since the controller will 
be self contained computer code communicating only with tue 
Senverters. Additionally the code will be faster since it 
eee NOt be necessary to Simulate the plant. 

mresentiy the realization of the code is not implevented. 
Mere software code needed for converters and the actual 
monitor needed for the realized hardware is not develored. 
Mes is not a major problem and should be relatively straight 


forward for & computer programmer with general tardware 
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experience. The actual control program will be identical and 
will only need to link into the new converter programs. [n 
addition the simulation display is not fully implemented. To 
make this project a powerful tool the simulation results 
should be presented in graphical form to the designer similar 
mothe plots at the end of this raper. If properly displayed, 
when the designer changes one part of the system all effects 
throughout the system should be dynamically evident. The 
advantages of a completely containe? hardware, software. and 
display system are obvious. 

The time required to compile the programs is tefinitely a 
function of the computer hardware. Most newer ZE&S based 
Single user systems can compile the code in abdout fifteen 
Minutes. &@8€ based systems, older ZEG based systems, and 
Multiuser systems that snare the central processor will 
momeere wp to thirty minutes to compile. skecall the compile 
time is not critical as mentioned earlier and the reélization 
eae will run in the order of milliseconas. As a final note 
the programs would run much faster if only one langvage and 
no preprocessor were needed. 

Another area that requires discussio2 is the theory 
mivolved in the design of the correction factor. TIevelopment 
Meeea «general prediction factor due to ‘delays caused by 
meme rsions and computer calculations is covered very briefly 
Mmeveis paper. The practical results of a predictor that is 


Meeo) With the example problem is shown in the graphs at the 





end of this paper; however, how the results were obtaine and 
what they really mean is not fully investigated. Hopefully, 
this will be pursued and discussed in a following paper. Thre 
predictor is more a control engineers problem and needs to te 
investigated extensively to determine the optimum predictor 
mor all possible delays. #urther discussion is provided in 


the example problem. 
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IV. FUTURE EXTENSIONS 
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After a careful investigation of the types of problems 
encountered with various hardware and software interfaces, it 
becomes very evident that a sixteen bit on larger 
microprocessor that offers high speed calculation irtegrated 
mreo a complete graphic and plotting machine will he a 
definite boost to the realization of future work, 
muemertonally a math processor integrated circuit that can 
operate concurrently with the central processor will Xe of 
invaluable importance. Once you place these bkardware items 
merece with good converter circuits into one unit it wiil no 
longer be necessary to simulate in one machine and realize in 
another. Once the simulation code works all that will be 
needed will be to specify executing the realization code on 
the same machine if connectors are available on the machine 
for communicating with the plant. Once the system overates 
correctly the nardware woula de reduced to a smaller special 
purpose system; however, the power of using a microvrocessor 
more tne controlier cannot be overlooked. If one vroviaes the 
Mepaoility of a full blown microcomputer in the sontrol loop 
even though all processes are not needed, the system could be 
sent to a distant changing environment and if the control 
equations needed to be changed due to unforseen circumstances 


it could be done remotely. By using the above metnod it would 
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be possible to reconfigure and realize a a2ew system 
dynamically which cannot be done by a physical one time 
hardware implementation. 

amonmer requirement 1s the software system needed for 
this type of development. The problems of integrating two 
completely different software packages is evident in this 
paper. An ideal software package would be a scientific 
relational database that has the capaoility to control the 
computer hardware and provide selt modifying code or 
automatic programming for it°s own operating system. This 
type of language would make the programming more efficient, 
ease the implementation, and provide a friendlier environment 
for future expansion. This type of softwere package is 
realizable under the new concurrent languages being developed 
such as ALA (product of the Lepartment of Tefence). ATA has 
not been proven to work in microcomputers since it 15 such an 
extensive language; nowever, several companies are showing it 
to be feasible in the newer computers such as tre €8@000 
Meecroprocessor from Motorola. The 68600 is a thirty tw. bit 
processor with a sixteen bit external bus. Motorola is 
presently designing a version of the 6&400 to aave a thirty 
two Oit external bus. 

It is recommended that future systems contain processors 
Such as the 68909 and math processors such as Intel's 8¢87 or 
Similar math processors from other companies such as Texas 


Instruments. Additionally high speed high precision analog to 
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digital ana digital to analog converters should be used. The 
type of converters used should contain their own multiple 
mopt Sample ana hold circuits. Tnis would insure capturing 
all inputs at the same time (versus time multiplexing) making 
tke Souverole realbizeations More practical for real time 


meeretions. 
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V. DESCRIPTION OF DATABASE ALGQRITHOM 


fee GENERAL DISCUSSICN (Refer to Appendix A) 

This section covers how the database turns the functions 
into PL/1-80 code. It sovers each suoroutine and there are 
hints on how the control designer uses the system. The 
meeeweSt way for the control aesigner to use this thesis is to 
m@mrerly read this materiel, go through the example ovroclem, 
and then proceed directly into the program. In geeneral a 
database is a program that contains related information. 
Normally tkese programs are used in business applications 
moore it is desired to keep track of an individuals name. 
address, or account information; however, cdatabases can be 
used to Store other int Onmat lon; Ibe this program 
implementation the database contains the function name, the 
type of function such as an analog function, names of other 
Memctions that are used by this function as inputs. the 
Seations for the functions, the relationship cetween 
memect ions, and if the function will be saved for later recall 
mao plotting. 

mieee prOsatam haS tae ability to develop new or recall old 
system models that have been previously develoved. Cnre the 
Peoeram is entered the model can te changed by addins or 
Memetite transfer functions, change relationships. or change 


meemed uations. édditionally the models initial conditi 
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Domechangead. It only the initial conditions are changed the 
program will compile much faster since each function will not 
need to be recompiled. After the system model is develcped 
the database proceeds to write PL/1 include files that are 
later included into standard PL/1 modules developed for this 
mmo ject. 

The features of the LEASE II program are the high level 
relational database language and the ability of the program 
Mmomewrite text files onto floppy disk. The disadvantazes 
include a limited math capability and it does not provide the 
ability to interface with the microcompute> hardware. A 
severe drawback in the program is that it cannot produce a 


memrcolon in any of the text files that 1t writes. 


B. LEVEL #1 (MAJOR SUBROUTINES AND PROCEDURES j 
Pee oNU.CME = initial menu selection 

Prmenidices dilevariadl2s and tne database. 

b. Enters a CASE statement to select a new proolem, 
old problem, save the present probiem, display ovpresert 
function names, change only the initial conditions of the 
meesent problem, or quit the program. 

foeeeoet eine S LOR and Lsk in Appendix A. 

meee OumNU.CME — second level menu selection 
a. Develops the control model database. 


db. Enters a CASE statement to edit, add, delete, 


meeall deleted functions, or change initial conditions. 





Ce Cates CUE wou. GMT . 
d. See Figures 16A and 165 in Appendix A. 
Oo. CONTROL.CME = database control program 
Ao ce omproceduresebOGVON.CMED, FROC.CMD, DATA.IML, 
and FORM.CME that develop the FL/1 include’ source code. 
b. Ends the database program and invokes the PL/1 
programs. 
eee LVON.GMI — monitor declaration program 


Sele Velous: UGCLMON. TXT for the MONITOR.PLI program. 


Demmodilns CNP GL. CML. 


(yn 


MEO NCO mom Gevelops AA.TAY, j|AD.TXT, OCD.TXT, ana 
meee. include functions for MONITOR.PLI. 
Seer eOcmonl >= Levelops PCLSUN.AXAX, INPUT .EXX, FUNT XXX, 
Mme OU. kAmM for e@ach function developed by the control 
Meme er to be used in the appropriate function HAA.PLI, 
meet, DD.PLI, or DA.PLi. 
Eee lale.Orl — data declaration procedure 
qe Leveloos BRE AIN¥O.TAT to be used is PARINT.2LI 
pe! ce Ce lt Ar contains Yaluveole format information 
pertaining to the LATA.FLT file. 
Cec cee omen ec ue — eaaditional data format f 
pees Oke Ci > PL/1 Command execution program 
Meevc vouonDGms omni, tlle 2O2".SUF used to control the 
Mecer Of FL/1 program compiling. 


Peeeiocwes me COmmands to reformat files conteinine a 


Meeeslesh into files with a semicolon. 
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Coeeeia liom UCLEUN-cL = function name setup for REFORM 

d. @@aills REFORM.CMD -— general format for FL/1 
@ompiling. 

Eee nOkiaenD — Short compile program 

a. Tevelovs initial conditions for TCLMON.TXT and 
MONITOR.PLI. 

b. Issues commands to control the compiling of only 
MONITOR.PLI by developing a submit file SFORM.CMI. 

Ceca | WcomormOnMecil = Short format for®PL/I compiling 


Similar to FORM.SUB called SFORM.SUBE. 


fee LEVEL #2 (PROGRAM DESCRIPTION) 

meccometune Majority of the code is in higher level 
languages, and is designed in modules, it is already in 
algorithmic form. [t is assumed tkat the individual needing 
to understand the code will be familiar with references 1 
through 4; therefore, an extensive discussion of the ccmputer 
code will not be provided. The example problem will cover all 
the features of the program. 

mne Dasic structure of tne database includes a memory 
file that contains the value of initial variables and a 
database if HALTS toate Contains the functions and their 
relationships to each other. ‘Tne entire database program is 
involved with manipulating and recalling information from the 
two above mentioned files. 

Mae initial memory file is shown in Figure ¢ and is 


Memeained in file INIT.MEM. This file is provided so that 
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Meytid: Variables Called will not produce errors in the 
database program. The main variables of concern are: 

ie tc = a@iscrete sample period, 

Pk ao aialog Sample period, 


See eoulMiT =~ end criteria for the digital simulation, 


ees oeMl rl = number of andlog samples for each digital 
sample. 
Bere LeLAY = delays due to calculation and conversion 


expressed as the number of analog counts after the sample 
period fT. 
Unseen system variables include: 

1. ACTION - users next menu selection, 


Zeon = NOldS database program name, 
ecu test Varlabie, 
Mover test variable, 


GS. ELEMENTS - data array position for a variable, 

faoouNe = index counter, 

epee UN > String of characters that is used to unigqnely 
define a function filetyve, 

EeeuniAl = number of datd array elements saved. 

mae Daslic structure for all functions is shown in Figure 
Seeeond iS contained in file INIT.DBF. The initial values of 
mee ddatabase are shown at the dottom of the Figure 3. The 


only items the control designer sees are listed under the 


NAME column of Figure 2. The user must enter values for the 
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Figure 2 


Used tone ini tialization 





STRUCBOGES FOR FILE: “INIT .IBF 
NUMBER OF RECORDS: COGG1 
DATE OF LAST UPDATE: €8/20/2¢ 
PRIMARY USE DATABASE 


FLT NAME TYPE WIDTH CEC 
OO1 FUNCNAME C 288 
O@2 FUNCTYPE C OBe 
003 NRS INPUTS C OO1 
O64 INPUT:A C 008 
O85 AR C 061 
O06 INPUT :B CG C08 
OO"? lal era 1B C O01 
088 INPUT :C C C08 
OO> CiE2) C O01 
010 Pee 20 CG OOS 
O11 DiPeL C 061 
G12 INPUT: & C GTS 
013 lopbie fale. C OO1 
014 FK N C08 GO4 
a) alee) Ge N BYE w04 
016 HQUATION CG G3 
017 SAVELATA L CO1 
G18 LATANAME C G85 


* Below is a list of the values used in the initialization. 
* These values correspond directly to the structure above. 
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variables under NAME of the correct type and size specified. 
The definitions for the variables are given below: 

1. FUNCNAME - function name used to identify this 
function initialized to “NEW FUNC. 

2. FUNCTYPE - describes this function’s type as: 

ae AA ~ analog, 
pee AYS= analog to digital, 
Clem L =mdael tal, 
dee DA ere tal to analog. 

Secwedmo. = TRissis the aumber of other functions 
mors function uses for inputs. AL(x) and CTA(x) registers are 
Met counted as functions. 

oe nis is a variable for name of function 
used as an input and is shown at tne bottom of Figure 3 45 
BONE . This valwe 1s referred to in the equations as IAP or 
Meyetor INPUT A PRES=SNT or INPUT A DSLAYED. 

a cee iS ls a Goede to tell tne computer the 
designer is talking about the “A input and he wants either a 

Mee for the PRESENT value or “2 for the TELAYSE value. 
matering a "4" will make beth the present and delayed values 
available but the program will run slower than if only the 
aqelayed value is needed. 

See nbss) = 18a5 variable is identical to 4 atnve 


Beeept this is the 8B input and has an initial name of TWO 


feesnewn at the bottom of Figure 3. Ps C.”)6CCUIINPUTEL, 6 | and 





INPUT:E are the same as above which allows for up to five 
maAputSs . 

eo “similis vdriable is the same as 5 above except 
mess tor the 8 input function. C1P2D, D1IP2ZD, and &81P2D all 
have the same function except they are used with their 
Mmeooective inputs. 

en om erase tunetion is the present value of the 
function named under FUNCNAME above. 

G. KD -=— this is the delayed or previous sample period 
value of the function named in FUNCNAME ‘is. £K delayed). 

moe cOUREION — this is the equation for the function 


named in FUNCNAME above. The equation takes three forms; 


Fo= cerm----- this is a function, 
AL(x) = ----- ValMewhor an analos to difzitval, 
DA(x) = ----— value for a digital to analog. 


meee ono bare — logical value ingicating if the user wants 
Memodave the value of this function for tabulation. 
Preun Anite meif it is decided to save the function then 
mat type tunction is its 
io—etne wna cCuion itself, 
DL xe = One of the eight analo,;- to digital, 
Peewee — One Or the Client digital to analre. 
It should be noted that functions are of three types. 
MemeecOmputer tinds a location for the functior or it is not 


hardware dependent then the function is described as FF. If 


the function is one of the other two types AI(«x) or La x) 





then they represent physical fixed registers within the 
computer. If the latter two types are used as arguments in an 
Saation they are not listed as inputs under NR:INPUTS for 
the function. Since they are physical registers they are 


called by their name (ie. AD(S)) and not by IAP or IAL. 


D. LOGICAL STRUCTURS 

The logical structure discussed will be a general outline 
See program flow. A complete discussion is covered in the 
example problem. Initially the control designer enters the 
program by responding [BASE MENU. to tne CP/M operating 
 —obem prompt. This will place the program under control of 
MENU.CMED which initializes the program and enters a UASE 
meaGement that provides a menu selection shown in Figure 15,4. 
At this point the designer usually starts a new problem Gr 
recalls an old problem that modifications will te mace toa. 
Notice that in order to save a problem the user must select 
Ppom this menu option. In general the next choice of the 
designer is to enter the modification phase either selecting 
memeechanse only the initial conditions of an old problem or 
selecting to make major changes. If major or new changes are 
mmesen the program transfers control to CMENU.CMD where the 
designer is allowed to change the function database by menu 
selections shown in Figure 16A. The choices are controlled 5y 
mie CASE statement in C™ENU.CMD. Most selections are single 
command entries followed by a RETURN from the xeybdoard. 


a) 


Some commands will operate without a RETURN especially if 





the user has exceeded a field in one of the database records. 
The fields are marked on the display with a =: at each end 
of the field. On some terminals the arrow keys may not always 
work for positioning the cursor in afield desired. In this 
case use CTRL E” (control key pressed at the same time as 
fae & key) to move the cursor up a field. Use CTRL X to 
move down a field and CTRL C to move down a complete 
mecord. 

Once the user is sure the database contains all correct 
functions and relationships he should return to tke MAIN 
PROGRAM if he wants to save the information otherwise any 
changes will te lost during compiling. The program could be 
improved by adding a question in tne CONTROL.CMD procedure 
asking the user if he wants to save the information o»oefore 
compiling. It takes more than fifteen minutes to compile the 
Meerram so a few extra minutes checking for correctness 
before compiling can be well worth the effort. 

Once the option to compile the program is taken, zontrol 
meses FO CONTROL.CML. This is the main control frocedure 
that causes the automatic ovrogramming oy calling all 
Mecessary procedures required to write the PL,1 sovrce code 
from the database. A complete discussion of this procedure is 
Contained in the example problem. After the PL/1 source voie 
is written, a COnUROl DEOgram 1S Written on how to compile 
the PL/1 code. Finally the database transfers control to the 


control program it just wrote, called FORM.SUB, by exiting 





mre Gdatad0asSe and submitting FORM.SUE to the CF/M operating 
system. At this point the PL/1 programs are compiled and 


executed. 
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VI. EESCRIPTION CF PL/1 ALGORITEM 


A. GENERAL DISCUSSION (Refer to Appendix B) 

This section describes the PL/1 programs and how they 
interact to provide simulation of a computer controller fer 
an analog plant. There are five basic procedures consistixg 
Of a monitor that insures all functions are calculated at 
least once during each sample period of the computer, an 
analog procedure for describing analog calculations. a 
digital to analog procedure and analog to digital proce.iure 
for describing conversions, anaa caigital orocedure for 
MeecriDing computer calculations. There is only one monitor 
but there are as many of the other procedures as needed to 
describe the entire system model for each function. As an 
example if there are two analog to digital equations in tne 
problem then there would be two analog to digital procedures, 
one for each equation. The example problem will demonstrate 
the function relationships. The basic procedures require 
macormation for each function and that is provided ty the 
control engineer through the database program in the form of 
Mmreuude files. Once the program is ready and compiled 
automatically by the database program a result file called 
ANSWER.COM is produced. The problem is run dy typine ANLWER 
Momeecone CP/M operating syStem Prompt. 4s the program is 


running intermediate results are displayed on the diszlay 
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Mievice for troubleshooting should the results not be as 
Eemected. After the program runs it produces a file DATA.FLT 
containing floating point data. This file can be translated 
into a text file by typing PRINT FILt where fILE is a name of 
peeile provided by the user that he wants the results pleced 
peel. O « 

PL/1 was chosen as the operating system language beceéuse 
of several features. First it is a common high level languaze 
mie is transportable to many of the present microcomputers. 
Memeeeintent of this project was to <eep the pragrams as 
machine independent as possible and using FL/1 source code is 
fueesood way to dothis. Pi/l also has very good math 
capabilities. Since the type of problems veing solved are 
engineering in nature it was necessary to have a high ievel 
math language. ha@ Vrlonvea lly, cicuroraiy Lor Fist of6 
communicate with the computer hardware is a definite asset if 
the programs are to be realized into hardware. sSeverai other 
languages were looked at such as PASCAL but PASCAL does not 
Beave the ability to talk to the hardware. The C # £leéngva¢ce 
MeeeatsSO §6l00Ked)6«6at «but it lacks the matk capability. 4s 
Stated in Chapter 1V¥.; hopefully, ADA will cverccme all the 
Mrraculties by combining all the needed features into one 
language. 

PL/1-S@ has some problems. There is a definite strobslem 
Memeoeconversion factors as was discussed in Chapter III. The 


Meee Droblems have to co with conversions tetween floating 





Pont and other forms. The original approach was to use only 
floating point values and therefore not need conversions but 
it was later discovered that several functions do conversions 
rexardless of starting and ending types. This is a pretler 
and cannot be overlooked by the designer if he has protiems 
Mmmm results. If @gquations requiring complex functions are 


used the designer will need to become familiar with FL/1. 
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Po uevehL #1 (MAJOR SUERCUTINES AND PHOCELJA: 


) 
1. MONITOR.PLI - main control loop program 
a. Causes iteration of ail functions. 
Deaelosubes ager ume tions are callea in proper order. 
Gee oeves Gata for later cisplay and enalysis. 
eee asic tOrmat sor all -anziloz functions 
ae Gemels all functions used as imouts. 
pemeecaveulavese= ana Treturas 1t s own prese1t or nast 
memes to all functions that call it. 
USieeMeSeeSVMincirohizgation of calculations with the 
eeeevem clocks through event counts. 
eeerrel wee Oasic “pernat “for’a@ll dnmalog te fdicival 
memetions and performs identical functions to those shown in 
meme D, and 2c above. 
leet o-wOdsic format for all aizital furctions and 
Meme orms identical functions to tzos® shown in 2 adove. 
Mel eee ol) = DaSioC Cormat for ell digitel to analog 


memeticns and performs identicai functions to trose above. 





G6. REFORM.PLI - corrects all improperly formed PL/1 
“include files by changing the backslash produced by thke 
BDASE II vrogram into semicolons. 

fe eR INT .PLI = wtility program Showing how to convert the 
floating point file DATA.FLT produced by the main program 
into a different printable text file for display on the 


system display. 


C. LEVEL #2 (PROGRAM LES CRIPTION) 

mence the Majority of the code is in higher level 
languages, and is designed in modules, it is already in 
algorithmic form. It is assumed that the individual needing 
Mom understand the code will be familiar with references 1 
through 4; therefore, an extensive discussion of the computer 
code will not be provided. The example problem will ccver all 
the features of the program. 

Poe OaSic Structure of the program uses a monitor that 
meaqmentially calis all functions within tke comvuter once 
marane €ach digital clock. The system clock is used as an 
meme COUNtTeEr to provide synchronization for each function. 
Tme functions are developed by the control engineer through 
the database program and each functional relebponsn 29 
Mareating which functions call each other must be ovoroverly 
provided. This is taken care of when the designer states 
memece functions are inputs to each other. A function can nave 


MmeseoutDput connected to any number of functions but a 
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function can only have five inputs not counting analog to 
digital, digital to analog, and it’s own delayed output. The 
monction definition allows for very complex multiple 
input/output relationships within any single equation. A 
method for correcting for computation and other delays is 
meovided in the correction factor CF in the monitor. This is 
a factor provided for the example problem; however, it is 
felt that this equation will work for any control system of 
the type requiring a steady state error signal of zero. The 
CF is provided but ¢an be e@liminated simply by not 
referencing it in any cf the equations. The programs provided 
are for simulation only but should be easily realized into 


control hardware since the design equations will not change. 


¢ 


Meee cOGICAL STRUCTURE 

moe program contains a number of proceaures that are 
basic howd  vroalems.. These were discussed anove. 
Additionally each problem has uniaue features that are 
provided by the control designer through the database 
programs. These features are coded into PL/1 ‘include files 
Semecead in dt compile time. Once the programs are coded all 
miormation as to which functions are connectec tocether will 
momeeorovided. The monitor is a PO WHILE LOOP that calls all 
memeoe CO digital and digital functions once each count of 
Mmeemeclock. At Gach count of the digital clock ( DCLOCK : any 
Beeered function will pverrorm it“s own equation evaluation if 


Meenas not alreedy been calculated during tnis clock and 





return the correct value to the invoking function. When the 
monitor calls a function it does not use the value returned. 
Paditionally there is a DO WHILE LOCP inside the main monitor 
moop used to call the analog functions the number of times 
Meecitied by the designer to occur vetween each digital 
sample period. This is the analog clock ( ACLOCK ). Also 
during the analog clock when the ACLOCK is equal to the DSLAY 
specified by the designer the digital to analog function is 
mevored, When a function is inveked, the first thing that it 
Mees 1s check to see if it bas any invuts. If there are any 
Meruts then the input functions are called to insure they are 
the correct value before the equation is calculated. ven 
though the monitor is a sequential process invoking each 
Mmimetlon, once a fumction is called it will invoxe all of the 
moots tunctions needed for equation Parc ley 10 0 - some 
functions could de called many times during each clock 
devending on how many other functions use it’s output. Once 
each TCLOCK the monitor outputs values that have been 
meaquested by the designer to save into a file. 

The first vrocecures called by the monitor are all 
fetal to analog functions located in a file written by the 
database program called “ad.txt . The functions take the form 
dummy ptr = SUNCTION_NAMS where FUNCTION NAME is the name of 
ae function and the file #FUNCTION_NAME contains Ee 
Meeessary information to represent the function. The example 


problem has a more detailed discussion. The first thing the 
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metal to anaiog functions do is call the analog functions 
that are used as inputs. If these analog functions have other 
analog inputs they call their inputs before calculating their 
output results. This process will continue until the analog 
functions have no inputs or their inputs come from a digital 
to analog register that will contain a value for a correct 
point in time. Similar processes occur for files ‘dd.txt, 
Memeixts, and aa.txt representing digital, digital to 
amadlog, and analog functions respectively. The crder of 
Syebine is controlled initially by the monitor and is 
extremely important for all processes to obtain the correct 
results. Direct coding of the processes would be much more 
efficient but it must be remembered that tnese processes must 
work no matter how- the designer initially defined the 
problem. This method cf problem solution has merits even for 
mee experienced programmer. tkven though the code is longer 
mean direct programming this type of problem will typically 
require one or two hundred files introducing the possibility 
more Syntax and other programming errors. 

cmt Our Main functions d@a.vli, ad.pli., dd.pli, and 
da.pii all have the Same general structure. Easically the 
mmemson is Ccalied by tne monitor or another function, the 
meme cion Calls any inputs, and then returns tke correct value 
of the equation at tne time the equation was invoked. This 
Smememeement places a process on the stack through the PL/1 


Saviins routine until all inputs are availacie eliminating 


f> 
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the need for developing a special process manager. Using this 
approach makes the programs easily adapted to concurrent 
Meocessors and languages such as AVA. The conversion routines 
meepli sand da-spli do not have a provision for using their 
delayed values since they would not be available in real 
hardware. This forces a true representation of the hardware 
and makes the designer provide intermediate holding functions 
if delayed values are needed. Additionally tne actual code 
for hardware will be easier to implement at a later date. The 
@enmer functions dd.pli and aa.pli have provision to return a 
previously calculated value if only the delayed value is 
desired. This allows the code to run mucn faster and provides 
more powerful equations within a single function. All four 
mpmecesses alklow for up to five function inputs and the 
meocedure knows how to call those inputs “by using the 
variable ‘nr_inputs to determine to which label in the 
memetion to jump. Basically all functions the designer 
implements are written into one of the above four procedures. 
AS an example if the designer specifies a digital function 
Called DIGIT, the database will write the include files and 
mien Compile the dd.pli procedure into a file TDLI.REL. Next 
Soe compiled file TPD.keL will be copien into DPIGIT.REL which 
Mere DOW contain all the necessary information for the 
mmection LIGIT. Additionally the necessary information such 
as the name LIGIT will te included into the monitor so that 


moeeenows tO Call the function once each clock period. 





Before any of the PL/1 programs can be compiled, all the 
text files written by DBASE II that contain a backslash must 
be converted to semicolons. This is accomplished by the 
procedure REFORM.PLI. This is a simple pro#ram that reads in 
a file, translates the character, and outputs to a new file. 
The files reading from and too are controlled ty the submit 
files FORM.SUB or SFORM.SUB written by the [BASE II program. 

The programs compile into a program called ANSWEn.COM 
under the control of the database program. ANSWER.COM is the 
file needed to run the problem. The designer only needs to 
place this file on a CP/M based system and type ANSWEn in 
response to the operating system prompt. This will run the 
program and send the intermediate results to the display and 
memgeeune tloating point data into a file called DATA.FLT. 

Once the program has been run it may be desired to nvlace 
mitemiLoating point data into a different form for display anda 
meemeysis., this project can be greatly improved oy vroviding 
miemresults on @ graphical disvlay device. Both tacsular and 
graphical results are supplied at the end of this paper and 
Mm@emereader can determine for himself which is the best 
fommod. One utility program provided is called PRINT Pll. 
7215 program reads in the floating point data ty using the 
Mmeeormatvion contained in the file DATAINFO.TXT as to row the 
date was stored and then sends the information into a text 
MeeeeeesSpecified by the user. This text file can then be 


printed if desired. 
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VII. EXAMPLE 


A. INTRODUCTION 

This part of the paper will make reference to all 
MemectS of the project and provides the easiest method for 
understanding the entire paper. [t should be noted that the 
descriptions as to how the data is entered and displayed is 
somewhat misrepresented since the entire project was designed 
for an interactive environment. The actual data entry at 
runtime is much more formal with data fields that allow only 
specific entry in specific areas of very specific type. This 
is not a constraint but an aid. 45 an example if the user is 
to replay if he wants to save the function results re can 
Oily answer with T, F, Y, orN for true, false. yes, or 
no. The answers can he in upper or lower case and if an 
improper input is made such as "7 the computer will BEEP and 
walt for a correct response. Additionally some information is 
Supplied in reverse yvyideo or highliehted. In the event that a 
color system is used such as the INTERCOLOR the information 
Meee be in different colors. The only way to fully appreciate 
fmemalds for the designer is to run the program. In order to 
meeriae a realistic test for this program the sample problem 
fmeermeinclitude a@ plant that is at the limit of stability and is 
mertially designed to be controlled by 2 bane bang 


Controller. The process will include a control loop that 
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feeds back on itself which provides a definite headache to 
computer proxrammers because ot the synchronization problems. 
fee PROELEM SETUP 

Assume we have a continuous me plant and it 1s desired 
to produce a minimum time response ripple free controller 
L(k) for a unit step input. keferring to Figure 4A below, 
L(k) is the computer transfer function, T is the sample rate 
Sree the input, LF is the delay from input sample at time T to 
output sample including all computation and conversion 
delays, and the clamp is included since the computer output 


is neld fixed between samples. Assuming that the clamp has a 


° 


tt) 
opal 


Figure 4A 


Computer and Plant Model 





very 


the digital system, 


M@eactical purposes the same as the input to the 


has deen 


Mmeritial assumption will 


conversions 


mais will 


error in reasoning will be easily corrected for later. 


should 


due to computer 
masrenificant compared to 
Ppamotuine K2(t) = 


amo letting the input to 


etiects of the clamp and computation delay, 


short conversion delay compared to all other delays 


the 


demonstrated by many authors [Ref. 5]. 


and calculations in the computer 


allow for simple design eauations yet 


be noted that most control systems ignore the 


Computation since the delays 


VHemthant soutpuc, K1i(t) = 


yal 


OMmpouy ol mine Clamp is fcr all 


clamp. [hrs 


The next 


be to ignore any delays due to 
(ie set [ = 


the 


delays 
are usually 
the sample rate. ) 

ay) Cute: (alae 
E2(t) the 


the plant be by lenoring 


aesereor a@iscre ve 


@rrerence equations can be developed d@scribding the plant at 


the sample instances 7. 
in 


G@rserete form: 


These eauations hav2 been 


the past by many authors [kef. 


ceveloped 


6] and are shown below in 


Xi(k+1) = X1(«) + T * F2(xk) (1) 
eae om See) + Ketek) = ((T * T) / 2] * B2(x; a) 
where; 

d/dt({Xi(t)] evaluates to 41(x+1) 

d/dt(x2(t)}] evaluates to X2(k+1) 

StvmesGesNext taSt 15 to control the plant descritet by 
mmemecquations in 1 and 2above according to the design 
reauirements for a minimum time response ripple free 





momo ller. ocoince this is a second order system it will take 
two time periods to obtain control. Assume the initial state 
vector equals zero (ie. X1(@) = X2(@) = @). Since we want the 
second time state with X1(2) = @ or no ripple and X2(2) = t, 
mee Input unit step, this leads to the set of equations below 


after some manipulation: 


i 


X1(2) Tos s2(9) + tT = a2(1) = @ Ge 
arco OT eT ye pen 2e( og) + (Tt * T)/2)*B2(1) = 1 (4) 
Solving equations 3 and 4 abdove simultaneously leads to- 


oa = 1/(t * 7) =) 


CaJ 
AV) ay) 
= 

tl 


iy es TE) (5) 

Pyeme@etrnitron ot the 2 transform, equations 5 and € and 
the fact that the desired output requires F2(k) = @ for k > 1 
leads to: 

aa 

mre e= 17(n * FT) -— [1/(T * T)}z ae, 

Now to find the Z transform of Hi(k) we start with the 
problem definition of zero initial conditions and a unit step 


mrput . 


At k = @ 


= i - @ = 1 

Meek = 1 

Bed) = U(i) - X2(1) 
mer = T)/2) * 629) 
ere 0) /2) 7 /(T * 7) 
=) 1 - 9.5 a Bee 





tJ 
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NO 
I) 
cS 
VW¥v- 
do 
i 
<i 
Ae) 
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nor K 2. i El(k) = @ 

The above set of equations leads to the Z transform: 

mz) = 1 + isa (8) 

eos) =) sez) 7 Bitz), the algorithm for the 
computer, we can now take the ratio of equations 7 and @ and 
Gross multiply then transform to find the difference 
equations needed. 

(1 + 0.52 "#82 (2) 5) (ees ee i ae "#21 (2) (3) 

ee meee eek —1)e= [1/(T * T))=[ei(k) - E1l(k=1)] (12) 

B2(k) = - @.5*82(k-1) + (1/(T * 7)]* (e1(k)-E1(k-1)] 0 111) 

Equation 11 is the one used for the computer simulation 
and shown in Figure 45 alongs with the rest of the model. 
mace it is not allowed to use Subscripted or arrays fer 
yariable names in this program implementation the following 


Changes have been meade in Figure 48 for variables in 


mmmatiors i, 2, and ili: 


mereriunction s2k. 


Mm e—1)) referenced as the delayed value of 22K or O in 


Mee function E2K. 


Peete — 9 Xi(%) or Xi(n) if nm different from & 
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PLANT 


Figure 4B 


Test Program Flock Tiagram 
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Xi(k-1) referenced as the delayed value of X1N or O in 
m@e function XIN and IAL in function XeN1l. 

Monee = A2(k) or X2(n) if n different from 

X2(k-1) referenced as the delayed value of X2(k) or “O° 
maecme function X2N1 and [AT in function Abe. 

UN = Y(t) 

mt => (T * i) / j 

fxeeeoae jmitial design the sample rate of the digital 
system T and the sample rate of the analog system at discrete 
points T1 were considered the same( ie. i = j). The only 
meason tor providing the possibility of calculating the 
mmelog values at different points 1s5S to see what happens to 
the plant between sample points of the digital system at fT. 
In most problems T and “i will be normalized to unity and 
“3° will represent the number of sample points calculated by 
M@eeepiant for @ach unit sample period T. 

Now to account for the combined delays in the discrete 
feeovem due to analog tc digital, digital to analog, and 
@emoutation delays. This function is called CF cr the 
Mmerrection Factor. The Cr in most systems is a major 
Semeraeration in initial equation development and usually 
Mmeeemoles Kalman Pilters or predictors. The approach taken 
for this thesis was that the average designer would not know 
Moeecould not estimate these delays due to system variables. 
mio means the simulations would not work as expected and 


Meeadiess to say neither would the realizations of the control 


©) 
To 





algorithms in the real hardware. The attempt was to force the 
designer to consider the delays in his design and possibly 
Provide a g00d predictor initially for his use. It was 
assumed that since this type control system requires an error 
mimemal that is zero when in control this would indicate that 
the signal is zero mean. Since the instantaneous input is 
random it was assumed the error signal was gaussian and 
provide a standard exponential predictor to the error signal 
to correct for delays. The predictor chosen was CF = “a 
where e is the natural log and £ is the combined delays. The 
CF is not a complete guess and further work will be vresented 
aed later date. The results of the CF to this particular 
example can be seen in the plots at the end of this paper. 
Further research needs to be done to see if this correction 
Moeees for truncation, roundoff, or even noise in the input 
Signal. The point is that some type of correction needs to be 
added and this one works very well in this case. [f the CF is 


not wanted the designer only needs to eliminate it from any 


of his equations. 


C. AUTOMATIC PEOGRAM (Refer to Appendix A and £) 
Mmeesaesjetern will initially enter the program hy entering 
"TBASE MENU and a RETURN to the CP/M overatine system 
prompt. This assumes the disk containing the [TRASE programs 
is in drive A and the disk Containing the PL/1 programs is 
in drive £”°. The system can be modified to run ona single 


Mon a5 was done in the original system but the disx must 
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contain at least one megabyte and should have directory space 
mor at least one hundred fifty files. Ey using the ovrogram 
setup requirements above the database is entered by menu 
selection “1 of Figure 15A to start a new program. This will 
imitialize the system and enter the menu selection of Figure 
164. Next if menu selection “5 is entered the initial 
conditions of Figure 2 can be changed to the ones shown in 
Figure 5A. This is done by the display menu CINIT of Figure 
memeeeNOLe that some parameters of Figure 5A are effected by 
returning to the main menu and Saving the program under’. the 
mame STRT, otners are developed as part of the program 
Meocess such as ACTION which changes each time a menu 
merection is made, and others such as CATA are developed by 
tne program itself counting the number of items that have 
been requested to save. AS can be seen from Figure 5A, 
ALIMIT was entered as 12 which means there will be 18 analog 
Pareulations for Gach digital calculation. Since the initial 
Meeren used tne same sample time for both tne digital and 
analog (isWweeol muse ve adjusted to €.1 as discussed 
Meeyviously since Tf is normalized to 1. After the initial 
Memmeitions ALIMIT, DBDLIMIT, DELAY, T, and T1 are entered. tke 
mmeeram returns back to CMENU of Figure 16A. 

Next assume menu selection “2 is selected and the 
Mimetblons are added. The structure will remain basically the 
Mame With only slight chanzes shown in Figure 58. The 


database for the problem will look something like that srown 





ACTION () oes 

PARAM CC SAGER T 

CN (Cc) 

FN (C) FUNCTION 
ELEMENTS (N) Q 
CLIMIT (N) 12 

Aon TT (N) 12 

Tal (N) 0.182 
FUN (C) ABCLEFGHIJKLMNOPQRSTUVWXYZ 
INI (N) 7 

Ny (N) 1.202 
PAY an 4 

DGEE (N) Q 


Figure SA 


Variables Used in STRT [Initialization 


RUG FOR FILE: aRT.EBF 
NUMBER OF RECORDS: COOSA 
COMO GAs te UPNAMS: 11/11/82 
PRIMARY USE DATABASS 


FLD NAME TYPS ea DTH DEC 
201 FUNCNAME C Za 
2B2 FUNCTYPE C Cz 
GOS Nie uit 5 C C1 
BO4 Hee Chie Bek C ZOE 
205 ewe aeb C CO1 
BOE INPUT :8 C Qe 
O07 IE 1aigra C O01 
J68 TNE? : C C 208 
AGS CLP aL C O01 
C1¢ JE RESARTAL 18 C C08 
Oil DB arab C CO1 
Che INPUT :E C ZGe 
8 20) EVP ey C CO1 
it FK N BOE B24 
ZAnD FAT N 008 CO4 
2aK6 EQUATION C ES8 
O17 Sey Ea Ts L ZO. 
g1s VTATANAME C CUS 
Grea ble howe se 


Semlcture fer SLR, 


a 
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fief ieure 6A — GT below. Keep in mind that the display looks 
mete different due to color, highlighting, or reverse video. 
meaoe onlymeone cecord at a time will be on the screen. 

The remainder of the menu programs will not be covered 
memce they are self explanatory. It should be noted that in 
order to save the problem the designer is working on, he must 
first return to the main menu and then save tne parameters. 
mee program will ask for the name of the file the user wants 
to save the program in. The example problem was saved in file 
STRT. If there are questions about the menus loox at 
meocedures MENU.CME and CMENU.CMD in Appendix A. After all 
functions are entered CORRECTLY, and the program is saved if 
desired, then enter "@ for menu selection in Figure 164A. 
mnis ae eee cause the program to enter the autoratic 
programming phase. (Take a break because it will take about 
fifteen minutes.) If tke vrogram fails to compile it will 
Meondably be due to improper naming of functions. Function 
mames Cannot be arrays. I[f the results are not as expected, 
ever though the vrogram compiles, then the equations are 
Mmeedably in error or there is a strange conversion in PL/1. 
If this happens then unfortunately references 2 through 4 
must be studied. 

wie compllinge bezins control passes to CONTAOL.SML. The 
Meemet «function called is DCLMCN.CMD which writes the monitor 
Mmeearations LCLMON.TXT for MONITOR.PLI and files 4A.TTT, 


Meee Leal i, and DALTIT sontainiag the procedure calls for 
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RECORD 86201 


FUNCNAME > UN 

PONGCLY Fs SAA: 

NR: INPUTS Oe 

GPNPOT : A SON 

41P2T7 LE 1 

PNeuTs B :TWO 

Bieee v a les 

Toi s C sTHAEE 
Creal Sale: 

CVewr: D ;HhOUR 
ere. a 

LNPOLs E oe iTV E 
E1iP2D a8 lhe 

FK > 1.0000: 
FKL > @.0000: 
EQUATION F=1 .2@000 
SAVELTATA re ee 
CATANAME an 


RECORD ©0202 


FUNCNAME ak 
FUNCTYPE 48 UO 

Nein Ful S dO) 
INPOT:A SONE 

A1P2r ae 
INPUT:8 ‘TWO 

Ea ee rae 

Nee C ot eee 
CiP2D ig 

INPUT: TC :FOUR 

Be 218 et 

[Neu | sFIVE 

ih jive 18 Seles 

FK >: €.6000: 
FKD ; 8.80000: 
SQUATION P=(AD(1)-AI(2))*cP 
Sa LA TA see 
DATANAME r 


Figure 6A 


Database For sxample Problem 
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RECORD @90993 
FUNCNAME ah ank 
FUNCTY PE see? 
NieseNe tw eo alle 
INPUT:A ree 
A1P2r Ss 
Nee Uy; B ev 
BiP2D ra lg 
Perl: C STARS 
Crea L ec 
Pewee: D HOUR 
ieee rsa 
INPUT: sFIVE 
E1P2C “1: 
FK 0.0000: 
FKD G.0@090: 
BQUATION eee pe) (Eat) =. 5*0 
SAVELDATA odes 
DATANAME oda) 


AKCCRD B0004 


FUNCNAME eee 
FUNCTY PSE oA 3 
NR:INPUTS Oe 

Ney SA :ONE 

A1lP2?t a 

CN PUT: 8 - YO 

Beee 2 ee 
INPUT:C sTERES 
Gasee fb sale 

Pew $ L :FOUR 

Bk 1 al, ple 

INPUT: & sFIVS 

ee ae rae 

RK : €.8000: 
FKT Q.2009: 
EQUATION sF=TA(1)*T1+0 
SG AT A rahe 
CLATANAME ¥ 


Figure 6B 


Latabase For Example Problem 
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RECORD 88005 


FUNCNAME >X2N1 
FUNCTYPE SAA: 

NR: INPUTS 31: 
INPUT:A SEN : 
AiP2D $2 
INPUT:B TWO 

B1P2Er ales 
INPUT:C sTHREE 
Cee L alee 
INPUT:. >FOUR 
D1P2D s1: 

INPUT :E sFIVE 
E1P2Er 25) 

FX : 2.89000: 
hap > O.0020: 
SQUATION eerie Ti = T1l/e+lArL*T1i+0 
SAVEDATA te 
DATANAME :F 


RECORD O06 


FUNCNAME Sard 
RUNCTYPE A: 

NR: INPUTS ae 

INPUT :A :UN 

EZ Ee ra oe 
INPUT:3 > TWO 

Biee) setms 

Pes C ibjabe ad 

ee 2 aes 

Eu [ S20Un 
Pee eons 

INPUT: & hes 
S120 cellent 

FX : 2.6000: 
FKD Q.000@: 
Bea TON sAT(1)=IAP 
SAVEITATA oe 
DATANAME orpen le) 


Pieevdse FOr Hxample Problem 





RECORD @8807 


FUNCNAME sAD2 
FUNCTYPE moll Be 

NRL NEUES one 
INPUT:A 2X2N1 
A1lP2r 2a 

Pueurs 8 sTWO 

Bore. oo 

ete Utes C et 6s 
Clea) a a 

Ne Or sc >FOUR 
Reece. ee 
INPUT: ore 
Bape somes 

FK >: €.@000;: 
FKD : 8.0000: 
EQUATION SADC} S140 
SAVETATA a he 
DATANAME ie )cs 


RECORD @8008 


FUNCNAMS Dav 
FUNCTYPE ae bs 

NR: INPUTS a 
INPUT:A AN 

A1LP2T mt 
INPUT:B ee NO 

BiP2ZD a ee 

BNEUT + C edi 2k; 
were L ie ae 

OSA BARS sFOUR 
LiP2r salts 

PNEUT se Ey E 

eee bs Oe 

FX oo. 0000: 
FKL  € 0000: 
NQUATION SO evi ae 
SAVELATA ees 
CATANAME 507s ie ee 


Figure €D 


Vatabase For Kxample Problem 
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Memon vores ee Figure 7. Next PROC.CMT is invoked for each 
Meincmncwnmmypes Of functions. Ihis orocedure writes the 
declarations for each function called DCLFYUN.XXX where XXX 
meands for which function number it is (ie. A, B, C etc. and 
mie last two XX indicate the function type (ie. 4A, AI, [T, 
Mey. PROC.CMD also writes the function input files INPUT.XXX, 
Mmedetion equation files FUNC.XXX, and the output files 
OUT.XXX that are used to indicate which element in the data 
array eack function represents. See Figures 8A through evel a 
moe next function invoked is DATA.CND. This procedure writes 
meee tlle VATAINFC.TXT that tells how many functions are 
maoeed in the floating point file [ATA."LT and the location 
of each element. See Figure 9. Next the program passes 
meeyrol to FOHM.CMD. This procedure writes the program that 
tells the computer what to do with all tne files that have 
Miob peen Written. Basically it rroduces a CP/M SUBMIT file 
meee tell the system to reform or corréct all files that have 
peeevacksiash into files containing semicolons, then it tells 
the system to compile each function and place it under it’s 
OWN unique function name. After compiling each function it 
memes them all together into one program called ANSWER.COM. 
see Figure 19. Finally after the database program writes thre 
mel)! program it passes control to the program it just wrote 
Ewe, creates the Simulation program ANSWER.COM = 4s BE baSee 


mentioned. 





Lo EGU. bx T */ 


UN ext entry(fixed(7)) returns(ptr), 
X1N ext entry(fixed(7)) returns(ptr), 
X2N1 ext entry(fixed(7)) returns(ptr), 
alee Srneetery nixed )) Seturns(ptr) , 
AT2 ext entry(fixed(7)) returns(ptr), 
DAG. ext entry(fixed(7)) returns(ptr), 
EX ext entry(fixed(7)) returns(ptr), 
Bek ext entry(fixed(7)) returns(ptr), 
ih 


TP iveat swetic ext init ( : 
Wimettoag. Static ext init( @.1¢00), 
Demi t fized Stetic init({ 1¢ 
homer ise float Static init( 192), 
Gqoray tloet static ini t( 4) 
iecar ext, 

2 da taneerc) dec iis, 6) 


VoONAeTXT = / 


dummy ptr = UN (code)\ 
dummy _ptr = XIN (code)\ 
dummy _ ptr = X2N1 (some) \ 

moe kL.LkT */ 
dummy_potr = AD1 Gecae)\ 
dummy ptr = ADZ (code)\ 

eo. eK? */ 
dummy ptr = EK ec cle 
dummy ptr = E2X (sode)\ 

roe A.TXT */ 
dummy ptr = DA1 (code)\ 
Figure 7 


POommuor LeClareattronsmand Functions 





/* DCLFUN.AAA */ 


UN sproc(code) returns(ptr)\ 

akc 1 

Mmeinputs fixed(7) static init(9), 

ONE EXiecmieryorised(,)) returns(ptr), 
TWO ext entry(fixed(7)) returns(ptr), 
THREE ext entry(fixed(7)) returns(ptr), 
FOUR ext entry(fixed(7)) returns(ptr), 
FIVE ext entry(fixed(7)) returns(ptr), 


mas) fixed(7) static init(i,1,1,1,1), 
Merioat Static init( 1.0000), 
O float static init( %.89023)\ 
dcl 
maat ext, 
2data( &) dec(i5,6)\ 


BeeeeouinP.TRT -— THIS IS PART OF ALL FUNCTIONS ANT WAS 
PoeGuatOne sat LeesiS CASE OF NUMERIC 
COMMEOL SYSTEMS. */ 
arc 1 
1 IA based(ia_ptr), 
EeltaAP float, 
Bent fl oeyt , 
1 IB based(ib_ptr), 
eeier float, 
eal t loat, 
mec based(ic ptr), 
eelor fioat, 
eercl float, 
1 ID based(id ptr), 
eee float, 
erry float, 
1 IE based(ie_ptr), 
Presi TiLoat, 
é IkT float; 
Pe SINPUT.AAA */ 


meput{1) = ONE \ 

Mmrut(2) = Two \ 

input(3) = TYREE x 

input(4) = FOUR \ 

input(5) = FIVE \ 

PoeeP UNC .AAA */ 

F=1.029@ 

/* OUT.AAA */ 

meca{ 1) = fic(? y\ 

Mme Skip list(“data( i) = °,F a 
miaure sh 


Hime tlon leclarations end Zaqauetions 





/* DCLFUN.BAA */ 


X1N sopoe code) returns (ptr) \ 
moc 1 
meigputs fixed(7) static init(@), 
ONE amt emtry(fixed(7)) returns(ptr), 
TWO ext entry(fixed(7)) returns(ptr), 
THREE ext entry(fixed(7)) returns(ptr), 
FOUR ext entry(fixed(7)) returns(ptr), 
FIVE ext entry(fixed(7)) returns(ptr), 
moo 5) fixed(7) static init(1,1,1,1,1), 
Merloat static init( @.9900), 
O float static init( @.0000)\ 
ac | 

1 dat ext, 


2@data( &) dec(i5,6)\ 


MeeeDCLINP.TXT - THI5 IS PART OF ALL FUNCTICNS ANT WAS 
Poet erOR eH THRSTS CASE OF NUMERIC 
COMEROG Stores. “/ 
ac 1 
1 IA based(ia_ptr), 
Ee lear ft Loa t, 
Bea float, 
1 IB based(ib_ptr), 
eeeer float, 
mmol tL Oat; 
imelG based(ic ptr), 
Eee Te@iait., 
ee icl float, 
ime based(id ptr), 
Pee float, 
eet float, 
1 IE based(ie ptr), 
eeiee float, 
melee float; 
fei NPUT.BAA */ 


input(1) = ONS \ 
moet (2) = TWO \ 
mrput(S) = THREE \ 
input(4) = FOUR x 
input(5) = FIVS \ 


/* FUNC.BAA */ 
F=DA(1)*T1+0 

fe OUT.BAA */ 
mara. <<) = ftc 
put skip list(’ 


Paure 25 


Function Teclarations and gquations 





Pe LCLEUN.CAA */ 


X2N1 sproc(code) returnos(ptr)\ 

ake | 

mm inputs fixed(7) static init(1), 

X1N ext entry(fixed(7)) returns(ptr 
TWO ext entry(fixed(7)) returns(ptr 
THRES ext entry(fixed(7) 


Returns GoUr 
returns(ptr 
seme erage Mes ay a 


FOUR ext entry(fixed(7) 
RIVE ext entry(fixed(7) 


es 

) 

returns(p tr) , 

) Ms 

) ae 
og’ 5) fixed(7) static init ( 
© 
g 


aa 
Pefloat static init( 2%.#2@@02), 
O float static init( @.@800)\ 
ie | 
ieaat ext, 


2 data( 8) dec(15,6)\ 


Pe eECLINP.TXT - THIS IS PART CF ALL FUNCTIONS ANLT WAS 
Tieeev Or Tae TEESIS CASE OF NUMERIC 
CONTROL SioLume , ~“/ 
dcl 
1 IA based(ia_ ptr), 
eeoearp float, 
POAT float, 
® 1B based{ib ptr), 
e LBP fl oat , 
reeset float, 
Mere based(ic ptr), 
eeerGe float, 
marc -tloat., 
ime based{id ptr), 
ee tfP float, 
meet float, 
1 IE based(ie ptr), 
PeeeeerP float, 
ees, float; 
PeSINPUT.CAA */ 


input(1) = XIN \ 
miput(2) = TWO \ 
Input(3S) = THREE x 
input(4) = FOUR XN 
maput(5) = FIVE \ 


/* FUNC.CAA */ 

mea 1) T1*T1/2+1AL*T1+0 

Pe OUT.CAA */ 

Mmea( 3) = ftc(F aN 

Memeo kip list(“data( 2) = °,? )\ 


Figure &C 


Function Teclarations and Ecuations 





y= OCLSUN.AAD */ 


ATI sproc(code) returns(ptr)\ 

ac | 

me inputs fixed(7) static init(1), 

UN ext entry(fixed(7)) returns(ptr), 
TWO ext entry(fixed(7)) returns(ptr), 
THREE ext entry(fixed(7 returns(ptr), 
FOUR ext entry(fixed(7 Fevuonowotr ), 
FIVS ext entry(fixed(7 Returns ( Dur), 


) 
)) 

Meo >) fixed(7) static init(1,1,1,1,1), 
45%) 
5%) 


1S 
merroat static init( %.@¢00), 
Be tloat static init( @.¢@00)\ 
@cl 
moat ext, 


Metta ( 8) dec(15,6)\ 


foe Geli NnrP. Txt — FTEIS IS PART OF ALL FUNCTIONS AND ‘WAS 
Pies tORTIHE THESIS CASE OF NUMERIC 
COMO l SroleMs . =/ 
arc 
ietA based(ia ptr), 
Pane float, 
mera float. 
fees Based(ib ptr), 
Pelee float, 
Matt float, 
Meio dased(ic ptr), 
Zeer float, 
BelGe float, 
1 ID based(id ptr), 
aloe float, 
mele’ float, 
1 IB based(ie ptr), 
eelee float, 
males float: 
PemeINPUT.AAD */ 


input(i) = UN \ 
maput(2) = TWO x 
iput(3) = THREE ‘ 
input(4) = FOUR x 
Poeue(5) = FIVE \ 
eeeeUNC .AAL */ 
AD(1)=IAP 


fm OUT.AAL */ 
ee, «CO 4) C= 06ftc(At(1))\ 
Bee Strip list(’data( 4) 


iI 
“ 

~« 

p 
sea 
~— 
t-— 
—— 
~—— 
~ 


Function Teclarations and Equations 


‘o> 
OD 





PeOOLCLFUN. BAD */ 


AI2 sproc(code) returns(ptr)\ 

gic | 

Meeinputs fixed(7) static init(1), 

X2N1 ext entry(fixed(7)) returns(ptr), 
TWO ext entry(fixed(7)) returns(ptr), 
THREE ext entry(fixed(7)) returns(ptr), 
FOUR ext entry(fixed(7)) returns(ptr), 
FIVS ext entry(fixed(7)) returns(ptr), 


Mees) fixed(7) static init(2,1,1,1.,1), 
Meee loat static init( 0.2000), 
Beeloat static init( 98.23980)\ 
ae | 
meaat ext, 
Pedata( &) dec(15,6)\ 


Bemeeou INP. TXT = THIS 1S PART OF ALL FUNCTIONS ANT WAS 
FIXED FOR THE THESIS CASE OF NUMEEIC 
COMBAOL SYSmave. */ 


fees’ based(ia ptr), 
< IAP float, 
e IAL float, 

m 12) pasea(ib ptr), 
poeee float, 
peer tloat, 

Meic based(ic ptr), 
meige fl oat, 
eeecl float, 

mee Dased(id ptr), 
Zee float, 
pbc t boat, 

1 IE based(ie ptr), 
eee tioat, 
eelet float; 

J~ INPUT.BAD */ 


input(1) = X2N1 \ 

moet (2) = TWO \ 

input(3) = THREE \ 

input(4) = FOUR \ 

input(&) = FIVE \ 

PoeeuNC. BAL */ 

AL(2)=IAL 

Paeour. BAL */ 

Mia 5) = ftc(AD(2))N 

Memeo list( data( 5) = °~,ar{(2))\ 


MMe eon beclaratilons and Hauations 


or, 
(0) 





Pep CueUN.ALL */ 


EK ‘proc(co@e) returns(ptr)\ 
ac 1 
mr_inputs fixed(7) static init(¢), 
ONE ext entry(fixed(7)) returns(ptr), 
TWO ext entry(fixed(7)) returns(ptr), 
THREE ext entry(fixed(7)) returns(ptr), 
FOUR ext entry(fixed(?7)) returns(ptr), 
FIVE ext entry(fixed(7)) returns(ptr), 
moms) fixed(7) static init(1,1,1,1,1), 
F float static init( 9%.@@09), 
QO float static init( @.€90@)\ 
lve 

1 dat ext, 


Bed@ata( 8) dec(15,6)\ 


PeeeDOLINP.TXT -— THIS I5 PART OF ALL FUNCTIONS ANT WAS 
Pier POnreink THESIS CASE OF NUMERIC 
CONTR@: SYSTSMS. */ 
acl 
1 IA based(ia_ ptr), 
£ IAP float, 
eelAn float, 
1 IB based(ib_ ptr), 
eal Tlodt, 
mer ticddat, 
Merc based(ic ptr), 
meeGe tloat, 
peeliclt float, 
teri based(id ptr), 
2 IDF float, 
eelLi tlodat, 
1 IE based(ie ptr), 
eelese float, 
Bele tloat: 
PoeeOnEUT.ALL */ 


input(1) = ONS ‘ 
moti c) = TWO x 
input(3) = THREE XS 
input(4) = FOUR . 
ieout(sS) = FIVE 


/* FUNC.ADL */ 

F=(AD(1)-AD(2))*CF 

/* OUT.ALL */ 

meat 6) = ftc(F aN 

Mer skip list(“data( 6) = °,F ee 


Figure &F 


function LTeclarations ana Equations 





y= LOLFUN.BLL */ 


B2xk sproc(code) returns(ptr)\ 
are 1 
Maminputs fixed(7) static init(1), 
EK ext entry(fixed(7)) returns(ptr), 
TWO @gt entry(fixed(7)) returns(ptr), 
THREE ext ie hl ot yy) ere uuinmsieper ) , 
FOUR emt emtmy(fixed(7)) returns(ptr), 
FIVE ax Entmy(fixed(7)) returns(ptr), 
mes) fixed(7) static init(1,1,1,1,1), 
F float static init( @.08090), 
Serloat static init{ @.9000)\ 
ae 1 

iedat ext, 


2data{ 8) dec(15,6)\ 


PeeUGLINP.TXT - THIS IS PART OF ALL FUNCTIONS AND WAS 
FIXED FOR THE THESIS CASE OF NUMERIC 
CONTHOE SYSTEMS. =/ 


1 IA based(ia_ ptr), 
~@ IAP float, 
men float, 

imeiB based(ib ptr), 
eter float, 
eel float, 

Were based(ic ptr), 
meses float, 
meetol tloat, 

ioe based(id ptr), 
emer float, 

* 10D fioat, 

1 IE based(ie ptr), 
eee float, 

2 IED floats; 
PoeeinrUT BLL */ 


Meput(i) = EK \ 
Maout(2) = TWO \ 
input(3) = THREE \ 
input(4) = FOUR x 
input(5) = FIYS \ 


Pema JNC.BLD */ 

=(2JAP-LAD)/(T*T)-.5*0 

feomeouT. ELL */ 

era 7) = ftc(F yk 

Meemeskip list(“data( 7) = ’,F eS 


Function Leclarations and squations 





Pe CCLPUN.ALTA */ 


[Al sproc(code) returns(ptr)\ 

acl 

Meminputs fixed(7) static init(1), 

me K Peleedtmwentixeod( 2) jeretucns( ptr), 
TWO ext entry(fixed(7)) returns(ptr), 
THREE ext entry(fixed(7)) returns(ptr), 
FOUR ext entry(fixed(7)) returns(ptr), 
FIVE ext entry(fixed(7)) eee 
fees) fixed(7) static init(1,1,1,1,1), 


Meetoat Static init( 94%.809@), 
Perloat static init( @.¢¢@a0)\ 
ac | 
Peydat ext, 
Bedata( 3) declis,6)\ 


PemeoGINP TXT - THIS IS PART OF ALL FUNCTIONS ANE Was 
FIXEL FOR THE THESIS cASE OF NUMEPIC 
GOveRCino rst EMS. ~/ 
rom 
feetA based(ia ptr), 
Mee float, 
meen float, 
mets based(ib ptr), 
eelee float, 
meow float, 
tio based(ic_ ptr), 
eel Ce float, 
Eich float, 
Meee based(id ptr), 
eeeeee float; 
Zeid) float, 
1 IE based(ie ptr), 
eeier float, 
< IED float; 
PeeeENEUT.ALTA */ 


mapm>tl) = nex \ 
Maput(2) = TWO \ 
iaout(5) = THRE ‘ 
input(4) = FOUR \ 
input(5) = FIV: ‘ 
eeuNC ALA */ 
TA(1)=IAP 


/* OUT.ADA */ 
Mee «86S) = ftc(TA(1))\ 
Mmeeskip list(“data( 98) = °,[TA(1))\ 


Figure &&H 


function Leclarations and Equations 


~} 
ay) 





7 Seen 7 


Pamete «ired= static init( 100), 

Roedeijde tered static Gnit( 38), 
npida fimed static init( 3), 
Dede cede sbatic init( 2), 
eriiweam Sta tlc Imit( 2), 
Miadcerixeds static init{ 1), 
name ( £) char(&) static in 


OLN oe 
“X2N1 a 
eh omer 
“AD ; ; 
“E2K a 
“DAL _ 


Figure ¢ 
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XSUE 

Remon VARS TXT AA TTT 
RiGOnM Alat AT ADIT? 
RoMOnM tia kT TL? 
RovOnMelAca x. LAs T TT 
PLI MONITOR 


Rent evi kas ei NevT 177 
RoVCkM LOGEPUN- AAA LCLEFUN TTT 
PIP FUNC. TXT=FUNC AAA 

fee Ory OuUle aaa OUT iit 

PLI AA 

PIP UN.REL=AA REL 


REFORM INFUT.PAA INFUT.TIT 
Ron UGERUN.BAA DCL2UN.TTT 
Piper UNC CIAT=BUNC .BaA 

RoeOnM OUT. BAA OUTS TTT 

PLI AA 

PIP XIN.H&L=AA.REL 


RereonM INPUT SCRA INPUT. TTT 
RoFORM DCLFIN.CAA DCLFUN.TTT 
PIP FUNC. TXT=FUNC.CAA 
ReeOnvwwour CAA CUT TTT 

PUI AA 

Pee A 2NEO REL AALG REL 


Peotone olen, INPUT TTT 
eons Loe ie wae LGLSUN. TTT 
Poe UNG TXT SeUNC AAT 

Ber Ory OULCAAL OUT. TTT 

rome eae k 

Pie ye bieek Sb H=n > RSL 


Roane wi Ne Ul. GAL INSU TPT 
ror Lous UNeeal LOLPUN TTT 
PIP SUNC .TXT=FUNC .BAD 

Poon OUL. Bat OUL.TTT 

PLI AD 

Pee cee ee. a 11 





Poder wee Den ee ee 
REFORM DCLFUN.ADD DCLFUN .TTT 
PIP FUNC.TXT=FUNC ATT 

oO rimouren rE OUL.T TT 

jo Bie 138 

PIP BK.REL=CL.AREL 


Powe eet BrP INEUr ToT 
DEeOnM [Ole bEL LPCLFUN.TTT 
PIP FUNC. TXT=FUNC.BDL 

Boor OUT.BDD CUL.TITT 

met LL 

PIF E2eK.REL=C0.REL 


meOony PUL ALA INPUT .TTT 
FomGnMere LY UN.s 0A LOLYUN. TTT 
Peer ONG. Tht —-P UNG ATA 

REFORM OUT.ALTA OUT.TTT 

PLI CA 

Pre eva he l=DA. RAL 


Pav er LN 
LINK PRINT 


MeN ANOWSR=MONITOR ,UN ,XIN ,X2N1 ,ADIL ,ADE ,DA1 , 5K 


Figure 1¢5 
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Pee reStT REGULTS ANT PROGRAM EXECUTION 

As mentioned previously, once the xzwgzau ANSWER.COM is 
developed, all that is necessary is to enter ANSWuF to the 
CP/M operating system promot. Figure 4B and 92 should te 
referenced during this discussion. Figure 11 is a partial 
printing of what will be on the display when the vrogram is 
mg@ewtineg. fhe printout is for the sample problem but with 
Pee CELAY set to zero and the LTLIMIT shortened. The first 
thing tc note is that the functions are not executed in «rder 
and some are executed several times during each clock. Tire @ 
Meetaily runs from PULOCK = @ to ACLOCK’ = 1. By associating 
the data element number with the orcer cf the function ‘ames 
in figure GS the user can determine which function is oeing 
Seeeculated or called. 34S an example the first element 
returning an answer in Figure 114 is element 4 which ecquates 
moat! in Figure S. The next elements are &, €, and 7 that 
Pequate to ALe, uk, and S2K respectively. An important ncte is 
Meemoesince there is no delay, function 3 (IA1) or the dizital 


Mmomeaneloe function is calculated or called at time @. This 


fee 0t be the case when discussed later on. The only 
meme vlon that has avalue is element i or JUN. There is n0 
memeeer changes in the date until DCLOCK = 1 Which is the 


mmemeeoime 45 ACLOCK = & for the second time. At this point 
Miemeanalos to dizital conversion has taken place for the 
MeemtLime, the computer does internal calculation, diital 


memanalioz is complete since tnere is no delay, and the analog 


ie 





meee starts to intesrate the signal. This is reaily the 
start of the sample since the computer is operating on tke 
Peaniiee conditions for the first ten counts of the analog 
mock. Sy the first sample period after tne initial values 
(PCLOCK = £4) the controller output changes tc -1 and the 
miamt Output is at the half way value of the input. Fy the 
feeond sample period after the initial conditions ({UCcLOCca = 
Sweeure plant output 15 equal to the forcine function input 
mmeetne controller output is @. It should ne noted tha’ in 
all programs the Zero point of the sample pericd always 
meers dat LCLOGS = @ so that the initial conditions ca: oe 
meepeaed. ithe adove example is a theoreticai non rea:izatle 
Semgueation of a control system with no calculation délays. 


The results obtained are identical HOrwcC bass iste: (TEXT oaon 


fis 


examples. 


taj 


igure 12 is the identical proolem discussed vitn the 
delay equal to zero but it shows what the display file loo«s 
Mmeeearter the utility program PRINT.FLI is run. The results 
Mmieerigure 12 are obviously more readable since the funeticn 
Meme and final value are not presented untii after all 
Saeculations are com ilete. 

meee o tS a Dartial printings cf the results Wren usins 
Meee ill example presented abcve with a forty percert delay. 
Semememeans there is four ACLOCK counts after the fYfcLccs 
Metore tne dizital to analogs values are available. The 


Mimemetant points to note are that the digital tre enalog 


ie 





maiues (DA! or function @) do no change until tke analog 
memgek ( ACLOCK ) 18 equal to 4 during each digital clock, the 
overshoot is only about 2.5 percent, and it takes five times 
memeeong to reach stability or the input signal. These results 
are much more representative of a real world problem ana the 
problems that would be encountered with a true herdware 


mopkementation. 





@ ICLOCK 
gee Gh G .CLOPLGE+OO 
deena 5) = 2. ASSCOVE+OP 
data( 6) = G . POSSOSE+LO 
eee. @  PLACVOE+2L 

@ ACLOCK 0 TIME 
dawwa ( 8) = Z .OLOSSBE+OR 
dranwant sd) = 1. @OOGABE+2O 
Jaa 2) = 2 .OSOGWGE+LO 
Gesical 9S) = v  OCOCLOE+UG 

lec Oem ye 
dae ( we) = 1.@930GGE+O0 
lever 2 Z.LUOCLCHOE+KL 
data( 2) = @.G0CCOBE+0O 
aval Bj) Z C2UOWOE TC? 

SAC LC Ge 5 SES 
data( i) = 1 @CCOWONL+9O 
deme s2)>= ¥  OOSSBOzE +2 
data( 2) = G  DOCCCOR+BO 
@ata( oS) c= @. @2OGCOE+OU 

TRPO LCC K 
data{ 1) = 1. @0CGAGE+39 
aati 4) = 1.290029 E+00 
Gana. oO) = D  OCOCOSE+CuU 
deste ( 5) = @ .OUSOCAE+0O 
dee Bie 1. G820002+C2 
dewalt 6) = 1. OCYOEO +62 
dana 7) = 1. SO@OCVE +49 

g¢ ACLOCK 1 TIM: 
data( 7) = 1. G2OSAPI+VL 
datva{ ~@) = 1.Q0O008092+u2 
detvau eet ) = 1. OGCZOVE+E@ 
data( 2) = 1. S20O0@-21 
data( 2) = 1.O000998-21 
data( 3) = 4,g9cS6S98-73 

tee LOOK PTs 
eho | ele 1. OLGAOGE+EE 
foe 2) = Zz. C2GCO8t-v1 
davat. 2) <= 2 .20C8O9RE-Y1 
data( 3) = ~CSCSGCH-Ge 

Ease oles Ve BESS 
Gasvaue 1) = 1.@2Z2009E+40 
data( 2) = % GLOOvOE-C1 
@aeeeat 2). = Z.uGCGOOeE-G1L 
aca co) = 4,40CSC68-92 

Pogeure “1A 
ANSweR.COM Runtime nesults 


N} 
C5 





SVAC LOG 
data( 1) = 1 
datal 2) = 1 
fatal 2) = iL 
eRe (Bi = y 
2 Lect 
data ( Le il 
asi bei 2h) 1 
fava 6) = Q 
diava( 9 5) = Gv 
12 bel 6) = Q 
puedo. — v 
aca, 7) = = 
0 ACLOC: 
fata.) = —1 
Meta ( ee) = 1 
gota | = 1 
faba). = Gg 
Gata 32.) .= Z 
data( 3) = U 
GC ACLOC 
Gatac i) = 1 
fata 2). = ih 
data( 2) = 1 
mattal 5 — 1 
Seve Loc 
daca i = 1 
data{ 4) = 1 
teas) = st 
data( 5) = 1 
Matar 6). = ? 
@ata( ¢) = 44 
elke. Ya Q 
CPC oGG 
mata on) >= C 
data{ @&) = 4) 
Gate, 1.) = 1 
data({ 2) = 1 
Saat 2) = it 
data( 3) = 1 
lea LOC 
eat = 1 
@ata( <<) = 1 
slelacul = 23) 3= il 
CSE So) 1 
Figu 
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ANSWER.COM Runtime Results 
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TIME=19 


UN = 1.@0800 
X1N = 1.20800 
XeN1 = C .500009 
AD1 = 1.92980 
ATe 7 0 .9BG0GO 
EK = 1.G@20OU0 
E2K = 1.828000 
tal = 1.090002 
TIME=2@ 
UN = 1.20080 
XIN = 2.500002 
een 1 = @.55506@ 
AD1 = 1.00000 
ATe = G .520808 
EK = C.520000 
nek = -1.060802 
Ie ak = -1 . 220000 
Tiwe=29 
UN = 1.288000 
X1N = C .O0GOUS 
KeN1 = 1 . GYCUas 
AT1 = 1.802008 
ill 2 = 0.520000 
i ¢ = C .500000 
wed = -1.0900000 
ral = -1.2d000¢ 
TIME=S2 
JN = 1 .G2O08O 
X1N = J. CdUUOO 
Xe2N1 = 1. O@02Oee 
ATI = 1.080000 
ADe = 1.200000 
aK = 2.200000 
ELK = G.UUGGGO 
DAl = Z.2O0O0O 
TIME=51 
UN = 1 .G3OUSD 
KIN = Z.20008O 
KeN1 = 1. O2UOHe 
slog = 1 .dOSA92 
END = 1.020002 
BE = Z OOLOLO 
nek = G .JCSOOP 
Tal = U .SO0000 


tere he 8 


PRevtePhi With Zero Lelay 


OD 
nm 





UN = 1.80000 
X1N = ZC .JQ0000 
ean L = G .O80d02 
AD1 = C . BLOG 
ATe2 ~ 0 .2O0CLO 
KK = QQ .200202 
Ke = 2 .200GO2 
ard = C .2OO200 
TIME=1 
UN = 1.8000 
X1N - ©. WSOOO 
X2N1 = 0. G0CUO 
AME = 2 .CL29000 
a = & .2O0CUO 
ok = 0 .DCUGGP 
nek = Z .JOOUd2 
IEF re = C . 002080 
TIME=9 
UN = 1.202000 
KIN = C.200000 
Xeni = 2? .C2OdH2 
euiel = GC .S80082 
fae = 0 .C2000¢ 
q = 2 .GCO0R2 
nex = J .0OGOBC 
Bue = C.2d0dOdL 
TIMB=1¢ 
UN = 1.08OU00 
KIN = C OLYGIO 
Xen1 = 0 .BCOUee 
nee = 1.882002 
ATz = C .232000 
BK = @.6703513 
hak = @.670819 
TAl = % .OSOaae 
TIME=11 
UN = 1. JSOdO 
X1N = Z .ce0OUe 
K2N1 = Z.OYCSd0 
Atl = 1.080002 
AT2 = 0.000000 
bk = 0.670319 
Eek = @.97031S 
LAl = 0 .O02000 


Paweure Von 


Fee NI orbit ewith Forty Percent [Lelay 


Hd 
Ou 





TIME=12 


UN = 1.002002. 
XIN = 2.00288 
KAN 1 = 2 .YUGOBGP 
AD1 = 1.@22820¢ 
Abe = 0 .2GCGAg 
HK = Meet are INS 
EOK = 0.670319 
ened: = CG. CUA0S¢ 
TIME=135 
UN = 1.020800 
ALN = C.d20800? 
KZN = 0G. COO0S0 
AD1 = 1.990002¢ 
ine = C .2600c¢ 
aX = C.07 2519 
Ben = 9.670519 
DAl - u.LEOEO . 
TIMB=14 
UN = 1.230009 
XIN = G.2E7CS1 
X2N1 = GO ViGSoe 1 
AT1 = 1.003000 
ee a EC .W22000 
x = @.670319 
wet = @.678051S 
pyrak = 9.67¢319 
TIMB=15 
Tal = 1.9¢9000 
X1N = O. Noa S 
Nal = 0.012406 
ATi = 1.2300 
ADE = J.O20C000 
BS = G.570319 
ee = Vig OC oS 
Dal = 2.670319 
TIME=4S 
TIN = 1.680¢00¢ 
KIN = Z.¢412¢2 
LOZ Na = i) Jor ele 
AP = 1.82200 
ADB2 = 0.322899 
wk = OY .2022 1¢ 
EK = -¢ .@SE455 
DAL = SU, Were ae' 


Sieure 12 


Poe liar Glew Norty Percent Lelay 





TIME=5¢ 


UN = - 1.080080 
KAN = 1.919908 
AQT1 = 1.@20000 
AD2 = 1 V1 GZ66 
WK = -2 .019903 
E2Kk = -@ .€14894 
LAl a = 026455 
TIME=51 
ON = 1.000200 
X1N = C26 cies 1 
X2N1 = NOPE ales, 
AL1 = 1.22000 
AD2 = to 216266 
eg = -0.0190E835 
E2K = -Z .0145¢4 
DA1 = Sey ocr OO 
TM = oe 
UN = 1.920200 
Yin = U2OLZo0c 
X2N1 = 1.024298 
Aw = 1.392000 
We = 120 TS260 
WK = -% 910963 
E2Kk = -% 0143894 
LAL = =n wge elo 
TIME=585 
UN = 1.80CG2C 
One = v -VO2eo0 
OZ al = 1.925%747 
AT1 = 1.000000 
AD2 = 1. OVe22e6 
a = -@.819¢33 
Bek = -&.0145¢4 
Ded = = eee so 
TIME=54 
UN = 1.02000@ 
XIN = 0.21176 
HeN1 - 1 O2e2c6 
eel = 1.008029 
Ae = 14 weecO 
UK = -9 .J1G9%3 
been = -¢ .014E94 
DAI = —{ .@l145s¢ 


tere 92 C 


aoe ith forty Percent Delay 


(1 
() 





NM Aa ss 


TN = 1.090000 
X1N = -% .@203518 
Nee NL = ie Coe: 
AT1 = 1.990000 
s\ ibe = oO H6Zeo 
EK = -€ 010503 
ELK = -@.8145S4 
eh = -%@ .014894 
TIME=50 
UN = 1.220900 
XiN = -¢ 801527 
X2N1 = Nn son tee. 
Ses! = 1.800000 
Age = LEO O1S) 
nk = =Vrnc oat 
BoK = -@.014894 
[Al SO -9 014894 
TIME=57 
UN = 1.000900 
X1N = -% .903297 
K2N1 = 1.924915 
AGE = 1.000880 
AD2 = ic keZco 
ed = =O loons 
R2«K = -0.0148¢4 
DAI1 = —7.G14894 
TIME=¢9 
UN = 1.@@000¢ 
X1N = -@ .40024c¢ 
X2N1 = wb .S9C45u 
ADL = 1.000002 
be = 1.0@03¢1 
sa84 = —LuoOUeec 
es ~ CV . 081s 
[al = Q.0012¢68 
TIME=100 
JN = 1. @Ygoee 
X1N = -2.00012¢ 
X2N1 = Z.90C431 
AD1 = 1.060202 
syle = @.S0045G 
BX = Q@.2003c8 
Oyen = € .JO00s56 
ical = 0.201268 


Pheure Lop 
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VIil. CONCLUSIONS 


aw @eeans  awew qwew GF a a 


Pee tlianyor "Point or this teesis was to cemonstrate the 
Meepeeoiity Of writings a computer osrogram that can ce used on 
fmmmcrocomputer by a novice orogrammer to aevelop coztrol 
meen models, simulation, and vtealization. This tyne ci 
System was not possible in the past due to system memory 
Mmearpeuati ons anc unacceptable software programs. "re approach 
tazen in tris project wes to use commercially availacle nigh 
meee roecramming languazes to reiuce the soding time. AS was 
Geesenssed jin tnis thesis, the laneuages are sot sufficien: in 
themselves and more tkan one lanzuaze was needed ts comvlete 
memmmeeronjes st. 1t 1S adpoerent that a language that kas tne 
features of a good relational database coupiel with extrer-ely 
meager matcsematical functions and the ahbility to communi-cats 
Memeemetce COmpuLter hardware is a necessity for this tyne 
eo o-m tO be a useful tool. 

Peememogetaace Of Tecosnizin= trat control systers can 2e 
Moeerec and manivulated in a datadasSe cannot se overliooke’. 
M—~EtGemr lc! 501s approas= ~t0 ‘ce .1m2t it is possitl2 to 


aeveloo a large databas= of control func 


7 4 
ere? 
Oo 
— 
Pa) 
+ 
a5 
(ap 
ct 
2 
o> 
res 
+ 
$2 
Sy 
it 


Meee eo ano sed like 4 library. Once tue database is formed 
Memeeould © manipulated oy artificial intellicence. Tre 
Memerver itself Ccouid te told to develon 4&4 corttrcl systen to 


emmermma 265, 7red Set or waveforms using the information sicree 





faerie east find a best fit solution. Missin2 parameters 
Could ve searched for, ocr vest fit solutions ceulad oessibly 
be provided as easily as using tne the present database 
systems to find phone numbers and names. 

brews: Swenificeant that source code was used for the 
[iMeeetmediate lanzuase. A special purpose program zoula nave 
merewricten treat proviged direct function calls to sresific 
mepgweare; nowever, by using high .evel lanzueke sour*e “ode 
mmemorozram will run on many different macnianes just oF using 
their resident compilers. Tris was demonstraved by 
@mwmspeonrting the »srograms to four different mashises with 
Meee cent operafine environments. ty using source code tke 
Mmool=m was Kept at a high level of adstraction. The vrozrams 
meemecomputer needed to write were similar to Enslish which 
memeec 861 25S5en the confusion while developing more shan ore 
mmaesec files and vrograms. 


me cClLmGiatebFeGane VEry apparent early in the orogran 


b4y 


Memeropment wes the need for a torrestion factcr to taxre into 
ferment the delays in tae system due to calculations and 
Memmpenrsions. fhis p~preblem is aadressed in meng moiels duit it 
is se@idom considered in most computer solutions since the 
Mmmevatlion time 15s usually insisnificant. It was unetermine 
Meee SOlvineg for all possibie delays was n2xt to irpossible 
especially when imrlemented since the nardware itself would 
Memmeerriabie. A correction factor is wrovided that works very 


Meme xitnh the example sroblem. Furtner inrvestiscation needs te 





Meme aone USing other functions to test the results; however, 
meeis felt that the factor will work with any control system 
designed to have a zero steady state error signal. The basic 
approach is to have the designer develop the model assuming 
there are no delays and then estimate the delays to be 
emeaed to the correction factor. If the factor does not work 
then try to develop one from other experience. I[f this fails 
then use benchside engineering, make estimates and watch what 
the results do. 

Appendix C contains several graphs of date covering the 
example problem. A variety of delays were introduced an: the 
Mmsu5+s with and without the correction factor are provided. 
The problem was initially designed assuming no delay and itnen 
Bomrected for with the correction factor. It should be noted 
that in all cases, even with delays up fe SS Fuecvew. 
acceptable resultS were produced within five times the 
miele! Sample period design. This is significant Since the 
uncorrected functions had severe problems even with the «mall 
delays. When using the correction factor the longer celays 
meeeected better ana it is felt that this is due tea 
fepputting a control signal closer tc tre time the intut was 
Peete to sample the result of the control. If a signal is 
meet and sampled immediately, the probability that ihe 
Meerpense tO that Signal is the one observed is higher tran if 
faeeeeone delay were encountered possibly resulting in severe 


Meo Crror. Sven thoueh the system was designed fora step 


Q) 
( 





input there is evidence that the correction factor works very 
well for a ramp input. This is also shown in Appendix C. 

As a final note the power of a graphical dynamic display 
for the designer is invaluable. Reading through many pages of 
data makes it very difficult to see where the significant 
Changes are being made. A real time plot will cut desiean time 


by an order of magnitude. 





CENT 
CMENU 
CONTROL 
DBASEAPP 
DBAS SMOD 
LBASESAT 
PARINIEX 
£ORMI 
MENU 
St 
XSUE 
EeOc 


UBF 
CML 
CMT 
OVR 
OVR 


OVR 


NUX 
CML 
Pik 
MEM 
COM 
CML 


APPENTIX A 


DPBASS II PROGRAMS 


CINIT 
CMENU 
PARAMS 
DBASESRO 
DBASEMSC 
TBASETTL 
SFORM 
INIT 
MENU 
SUB. TL 
SHORT 
CATA 


Da 
al 
CBE 
OVE 
OVE 
OVR 
SUB 
DBF 
(Noe 
COM 
CML 
CMD 


CINIT 
CMENU 
PAaTAl 
Peon sie OL 
TPFASEMSG 
CBASEUPL 
CCLFUN 
INIT 
MENU 

Ze 
CCLMON 
FORM 


ig is 
ite 
CML 
OVA 
COM 
OVT 
CMD 
MEM 
LPR 
COM 
CML 
CMD 


CINIT 
CMENU 
TBASE 
PBASEMA I 
PEA ones 
DBINSTL 
FORM 


>; MENY 


Splecgat 
ZIPIN 
MONICL 
REFCRM 


Programs Required for [atabase Operation 


DER 
BER 
LOIS 
COVE 
OVE 
GON 
SUB 
CMD 
List 
CCM 
Ce 
com> 





* MENU.CME 
SET TALK OFF 


RESTORE FROM INIT * See Figure 2 
Sok INIT * See Figure 3 
SET FORMAT TO MENU * See Figure 15A and 15B 
CLO WHILE T 
READ 
ro CASE 


4 


CASE ACTION = “@ 


ACCEPT “DO You REALLY *WANT ToeOQUIT? Y/N 


miacy = —% 
QUIT 

ENDIF 

CASE SCTION = “1 
RESTORE FROM INIT 
USE INIT 
COPY TO PARAMS 
USE PARAMS 
DO CMENU 

ONS er ACTION = 2 


i 


ACCEPT “ENTER NAME OF OLDE PROGRAM FILE 


RESTORS FROM SPARAM 
USE SPARAM 
COPY TO PARAMS 
USS PARAMS 
LO CMENY 
ask ACTION = “3 


4 


ACCEPT “ENTER NAME OF FILE TO SAVE 


SAVE TO SPARAM 
COPY TO SPARAM 

Pol ACTION = “4° 
LISP ALL FUNCNAME 
ACCEPT “PRESS RETURN TO CONTINUE 

mise ACTION = “5 
ACCEPT “ENTER NAME OF OLE PROGR&M 
RESTORE FROM SPARAM 
USE SPARAM 
COPY TO PARAMS 
USE PARAMS 
LO SHORT 

ENICASE 
ENLIO 


¢@) 
tO 


i 


A 


TO CN 


Ora 


PAE MA davon TN 


Crea Ay 





WOO DWDDW WADA DWM % 


Secile MENU *** 


ENTER 
ENTER 
ENTER 
ENTER 
ENTER 
ENTER 


OMe OOM PR & 


sae) 
TO 
TO 
TO 
TO 
TO 


QUIT 

START A NEW PROBLEM 

RECALL AN OLT PROLL&™ 

SAVE PRESENT PARAMATES 
DISPLAY PRESENT FUNCTIONS 
CHANGE [INITIAL CONDITIONS ONLY 


EG. Se) 


eC 
eye 
S alte 
Ee eC 
erat 
1 Suged *7 
one O 
14,52 
er 25 
wee o 


SAT 
Dink 
SAY 
DA Y 


( This will allow the pvrogran 
to compile much faster ) 


ENTER ACTION #ACTION 


ear e ben 

MENU 
iimn@ TO QUIT 
ENTER 1 TO START A NSW PROBLEM _ 
"ENTER 2 TO RECALL AN OLE PROFLEM” ~ 
"ENTER 3 TO SAVE PRESENT PARAMATERS” 
LENTER 4 TO DISPLAY PRESENT FUNCTIONS ; 
ENTER 5 TO CHANGE INITIAL CONDITIONS ONLY 


_( This will allow the program’ 
to compile much faster ) 
ENTER ACTION 

ACTION 


Figure 152 


MENU Formet 


(O 


CA 





* CME 

ont eek 

eo Ow 
KA 
DC 


ENT 
ENITO 


MU. OMT 
ORMAT TO CMENU * See figure 16A and 163 
Pos . 
1% 
CASE 
CASE ACTION = “@ 
PAG 
INDEX ON FUNCTYP®S TO PARINITEX 
US& PARAMS INUTEX PARINUEX 
DO CONTROL 
PASteAGTIONt= “1 
ERASE ; 
TISP ALL FUNCNAME 
ACCEPT “°ENTER RECORD NUMEER TO EIIT ° TO CN 
ELIT &OGN 
ers ACTION = “2 


7 


r 


¢ 


APPEND FROM INIT * See figure 3 
pet 4 

Gist ACTION = ~3~ 
ERASE 


LISP ALL FUNCNAME 
BECHPT “ENTHR RECORD NUMBER TO DELETE ° TC CN 
Dior e HeECORL SCN 

Mis ACTION®= “4° 
ERASE 
CISP ALL FUNCNAME 
ogee! SaNTER RECORD NUMBSEx TO RECALL ° TO CN 
RECALL RECORI SCN 

oo ACTIONS = “5° 
Si. FORMAT TO CINIT meSece tiene 17! “and 178 
REAL 
Soler Ongar LO CMENU * See figure LEA and 168 
ACCEPT “’PRESS RETURN TO CONTINUE ° TO CN 

Mon eACTION = °6 
SET FORMAT TO MENU 
RETURN 

GAS 3 


tt 


See figure 1SA and 15S 


xTURN 





aone= file CMENU *** 


ENTER @ TO COMPILE A PROGRAM 

ENTER 1 TO ELIT TRANSFER FUNCTIONS 

ENTER 2 TO ATT A TRANSFER FUNCTION 

ENTER 3 TO DELETS A TRANSFER FUNCTION 

ENTER 4 TO RECALL A [TELETED TRANSFER FUNCTION 
BNGER S TO CHANGE INITIAL CONDITIONS 

ENTER 6 TO RETURN TO THE MAIN PROGRAM 


ENTER ACTION #ACTION 


Figure 16A 


* CMENU.FMT 

es ,15 SAY SNTEE 
@ 4,15 SAY “ENTER 
@ 5,15 SAY ENTER 
@ 6,15 SAY “ENTER 
@ 7,15 SAY "ENTER 
@ 6,15 SAY ENTER 
@ 9,15 SAY “ENTER 
fee, 22 SAY ENTER 
@ 12,35 GET ACTION 


TO COMPILE A PROGRAM’ 

TO ELIT TRANSFER FUNCTIONS | 

TO ADD A TRANSFER FUNCTION 
NO Cpe Lor se A Ta Awe ER FUNCTION 
Aedes A Oe n TRANSEER FUNCTION 
TOLCHANGS INITIAL GONDITIONS © 
TO RETURN TO THE MAIN PROGRAM. 
CP TON 


me OOP CAM FR © 
3 
ee 


Figure 168 


CMENU Format 


(OD 
cn 





oe File CINIT *** 


DLIMIT #DLIMIT (This is the number of times the computer 
wall do calculations. ie. The number of 
Sample periods.) 
ALIMIT #ALIMIT (This is the number of times each analog 
Calculation will be done during each sample 
period. Bffects how accurate the analog plant 
is modeled. ) 
Poe 8#DELAY (This is the number of analog calculations done 
before the digital to analog is availible. | 
it #T (This is the sample rate of the digital system. 
mem Must provide it 4f you are eoing to use it 
in your equations.) 
at #T1 (This is the sample rate of the analog system. 
It is usually the same as T above during desien: 
mowever., for simulation it usually works out to 
be T divided by ALIMIT since you usually want a 
more accurate plot of the analog plant to see 
what happens between the digital sample freriods. 
MOUMMUSteProvige ld if you use it in your eoua- 


tions.)(This is effectively double rate sampling. 


Figure 17A 
CaN EL 


~ a” 





DOWDDD DDD ODD D WOH DW DDD DDDDWDADDAD ED *K 


“DLIMIT™ 

Diy y 
(This is the number of times the computer 
“will do calculations. je. The number of 
“sample. pes ods...) 

ie ie 

eLIMEs 
e(Geeis is the number of times each analog 
“calculation will be done during each sample 
“period. Effects how accurate the analog plant 
bis Mowe ledr he 

"DELAY 

DELAY 
"(This is the number of analog calculations done 

“pefore the digital to analog is availible.) 


T 
“(This is the sample rate of the digital system. 
You must provide it iyo areeescine to wee it 
“in your equations.) 

te 


et: 

“(This is the sample rate of the analog system. 
PemecoeUoiaiiy —Thessane as Tiabove during desizn, 
“however, for simulation it usually works out to. 

“pe T divided by ALIMIT since you usually want a 

Puogeeaceurate plot of tne analog plant to S52e 
what happens between the digital sample periods. 
“You mus t provide iol rer ou Use -1t> ine your egies 


“tions.)(This is effectively double rate sampling.) 


Prenure 178 


CLINT. forma t 


a 





= CONTRO Rec Mil 

[TO LCLMON 

SOs Oo OneNtE 

LOCATE FOR FUNCTYPE = 
DO PROC 

OTORE @ TO INDY 
LOCATE FOR FUNCTYPE 
YO PROC 

OTORE @ TO IND 
LOCATE FOR FUNCTYPE = 
YO PROC 

STORE 9 TO INI 

LOCATS FOR FUNCTYPS = 
LO PROC 

TO DATA 

TO FORM 

QUIT TO “SUBMIT1 FORM’ 
SET TALK ON 

RETURN 


<@) 
‘ae 





= DCLMON.CMD 


SET ALTERNATE TO B:ICLMON * 


SET ALTERNATE ON 

LOCATE 

[TO WHILE .NOT. EOF 
aad “+FUNCNAME+ ” 
CONTINUE 

ENDLO 


COUNT HOM SevEDATA TO DATAL 


iedat ext, 


IN VN VOY 
. 8 SN » 


SET ALTERNATE OFF 
STORE “AA TO T:NAME 
DO MONDCL 

STORE “AD TO T:NAME 
TO MONDCL 

Sone LL TO T:NAME 
LO MCNDCL 

Sees LA TO T:NAME 
LO MONICL 

RETURN 


weeeNCCL.CML 


SHT ALTERNATE TO B:&T: 


SET ALTERNATE ON 
meer e FOR FUNCTYPS = 
Ir SOF THEN 
is 
KLSE 
LO wHILE .NOT. EOF 
2 @urmy ptr = 
CONTINUE 
ENDDO 
ENT It 
SET ALTERNATE OFF 
RETURN 


NAMs 


T:NAME 


oe 


i/ @ 


“+¢funcname+ 


(0 
ce) 


te 
° 


pedlatat totr(DATAL,5)+ ) dec(15,6)’ 


not needed 


“Neodeyes 


not needed if single disk 
ext entry(fixed(7)) returns(ptr), ” 


Here oat Shatic Ext init( +str(T,8,4)+ ), 
Mimeitoatl Siuatic ext jnit({ ~ +str(T1, 8,4)+") 
pee fired static init( + ST rtTEIMIt 4) 
Meuwitierioabe Static init( + str(ALIMIT,4) 
delay float static Viet teste (LeLAY,4)+~ 


Hy,’ 
+*), 
) 


9 


i Sen eee aS ik 





+ PROC .CML 
[TO WHILE .NOT. EOF 
SEFORE IND+1 TO INI 
STORE “CCLFUN.°+3(FUN,IND,1)+FUNCTYPE TO T:NAME 
SET ALTERNATS TO B:&T:NAMB * EK: not needed if single disk 
SET ALTERNATE ON 
? FUNCNAME+’sproc(code) returns(ptr)\" 
“ael © 
‘nr_inputs fixed(7) static init( +NR:INPUTS+”), 
INPUT:A+’ ext entry(fixed(7)) returns(ptr) 
INPUT: B+” ext entry(fixed(7)) returns(ptr) 
INPUT:C+° ext entry(fixed(7)) returns(ptr) 
INPUT: D+° ext entry(fixed(7)) returns(ptr) 
INPUT: E+” ext ry suatio init os aaa ) 
ve b 


¢ 


rd 
4 4 
¢ 4 


YA 


is 
is 
is 
i 
if 
? 
ik 
ik react; 


meed{(5) fixed(7) static ini S22) oa “+h1 
°, +C1P20+", “+D1P2I+°,  +E1P20+’ - 

Sroms STR(FK 96,4) TO TEMPOUT 

mee hep ioat Static init(°+TEMPOUT+’),’ 

Some STR(FKL,S,4) TO TEMPOUT 

MmeOurel oat Static init(“+TeEMPOUT+’)\" 

oe dc] ” 

mee i dat ext, 

» * 2 data( ’ Misras. Sy odie 1536) \- 

SET ALTERNATE OFF 

SPORE  WNPUT,°+S(FUN,IND,1)+FUNCTYPE TO T:NAME 

Se ALTERNATE TO B:&T:3:NAMB = Bs not needed if single disk 

Ser ALTERNATE ON 


? “input(1) = “+INPUT:A+°\~ 

? “input(2) = “+INPUT:2+°\" 

fee ineut(ds) = ~+INPUT:¢c+’\~ 

Meerut (+) = “+INPUT:I+’\’ 

fe input(5) = “+INPUT:E+’\~ 
E 


SET ALTERNATE OFF 

STORE “FUNC. °+$S(FUN,INDL,1)+FUNCTYPE TO T:NAME 

Set ALTSRNATS TO B: ST: NAME S2eiet meeded 1: single disk 

SET ALTERNATE ON 

? EQUATION 

Sor ALTERNATE OFF 

STORE “OUT.°+9(FUN,IND,1)+FUNCTYPE TO T:NAME 

oe ALTERNATS TO Bs&T:NAME “ be not needed if single disk 

SET ALTERNATE ON 

im SAVELATA 
STORE ELEMENTS+1 TO ELEMENTS 
meoceica( +oMmR (ELEMENTS ,3)+°) = ftc(’ +EATANAMES 7) \” 
meomtesmipoeiist( data(|/+STR(BLEMENTS ,3)+1) = °,] 
2? TATANAME+ “)\7 

ELSE 
? 

ENDIF 

Sef ALTERNATE OFF 

CONTINUE 

ENLLO 
RETURN 


100 





* DATA.CMD 


a, 


SET ALTERNATE TO ESLATAINFO * EB: not needed if single disk 


SET ALTERNATE ON 


? *timel fine Static imi t( PORE @iL Mi oAbPMiT,4)+ ),- 
Mecmt data fixed static init ( "+STR(TATAL, rot) 

MouNT FOR SAVEDATAP AND. (F UNCTYES =the TO) IND 

? “nr aa fixed static init(“+STR(IND, aie yee 

COUNT FOR SAVEDDATA .ANT. (FUNCTYPE = 2) TO INT 
feeaoead fimed static init(“+STR(IND, 2) hoe 

COUNT FOR SAVETATA .ANL. (FUNGTYPE = “) TO IND 


? ‘nr dd fixed static Capt Ar 2)" ye 
COUNT FOR SAVEDATA .AND. (FUNCTYFPE = TA“) TO INE 
Weeemr da fixed static init( “+STR(INI,2)+’),’ 
? ‘name (* +STR(DATA1, 3)+’) char(&) static ae 
STORE “AA’ TO TINAME 
Poe DATA 
LOCAT® FOR FUNCTYPE = ‘AA’ 
IF .NOT. EOF 

2? 
ENDIF 
STORE “AL” TO T:NAME 
BO WATA1 
moeeis FOR FUNCTYPE = “AL’ 
fee NOT. EOF 

— 
ENDIF 
more “If° TO ™:NAME 
POM CATA 1 
MechTE FOR FUNCTYPE = ‘DD’ 
Pieeor., EOF 


STORE “LA” TO T:NAME 


SET ALTEXNATS OFF 
RETURN 


* DATAL.CM 
Moet FOR SAVELATA .ANLZ. (FUNCTYPE = T:NAME) 
DO WHILE .NOT. EOF 

2 [{°]+FUNCNAME+[7] 

CONTINUE 

Mm .NOT. sor 

a, 

ENDIF 
ENDDO 
RETURN 


i Geal 





* FORM.CML 
SET ALTERNATE TO FORM.SUB 
SET ALTERNATE ON 


a Be * B: not needed if single disk 
fm xXSUB’ 

ee eeMOoRnM AA .TXT AA.TTT ” 
fae ORM ALD.TXT AL.TTT” 
MORE EORM DD. TXT DD.TTT’ 
mee REFORM EA.TXT LA.TTT’ 
? “PLI MONITOR’ 

STORE “AA” TO T:NAME 

CO LTCLYUN 

Sogn AC TO T:NAME 

DO DCLFUN 

eons CL TO T:NAME 

LO CCLFUN 

STOR® DA TO T:NAME 

LO ECLFUN 


i Ll PRINT’ 
fm GINK PRINT’ 
? “LINK ANSWER=MONITOR’ 
LOCATE 
DO WHILE .NOT. EOF 
2? °, +TRIM(FUNCNAME) 
CONTINUE 
ENDDO 
? 
SET ALTERNATE OFF 
RETURN 


| DeeruN.CML 
LOCATE FOR FUNCTYPE = T:NAME 
STORE @ TO INT 
moe metLe .NOT. EOF 
STORE IND+1 TO IND 
mrOne S(FUN,INL,1)+FUNCTYPE TO EXT: NAME 
DO REFORM 
CONTINUE 
ENTLO 
RETURN 


* REFORM.CMD 

? “REFORM INPUT. °+EXT:NAME+” INPUT.TTT ” 

2? “REFORM LCLFUN. “+EXT:NAME+° CCLFUN.TTT” 
2? “PIP FUNC .TXT=FUNC. “+EXT: NAME 

2? “REFORM OUT. “+EXT:NAME+” OUT.TTT” 

? “PLI “+FUNCTYPE 

eel + TRIM(FUNCNAME)+° .RoL=°+FUNCTY PB+”. 
RETURN 


MaKe 


28, 
tc 





eoHORT .CMD 


SET FORMAT TO CINIT * See figure 17A and i7B 
READ 
Sor TALK OFF 
SET ALTERNATE TO BsDCLMON *B: not needed if using a single disk 
SET ALTERNATE ON 
moon | & 
beomwnliLe .NOT. EQF 
a “+FUNCNAME+” ext entry(fixed(7)) returns(ptr), 
CONTINUE 
ENTLO 
COUNT FOR SAVELDATA TO LATA1 
T float static ext init(’+str(T,8,4)+ ern 
felt loet Static ext init” +str(Tt, Q@,4)+") 
DriMitetixed static init(° + str (LLIMIT, 4) 
ALIMIT float static init(” SeourVALIMET ,4) 
delay float static init (“+ Sr epELAY, 4)+ ‘ 
1 dat ext, 
? 2 data(“+str(LATA1,3)+") Nec S56) 
So eeaALLTeSRNATE OFF 
SET ALTERNATE TO SFORM.SUB 
SET ALTERNATE ON 
? XSUBS 


any 7 
ei) af 
) 


’ 


VINIVVIN YO 


.» N .» NN . 8 


ii * BRB; not needed if using a single disx 
? PL MONITOR“ 
Zeit PRINT’ 
? “LINK PRINT’ 
2? “LINK ANSWER=MONITOR’ 
LOgaT E 
TO WHILE .NOT. EOF 
2? “, +TRIM( FUNCNAME) 
CONTINUE 
ENLLO 
; 
SET ALTERNATE OFF 
Gemt TO “SUBMIT1 SFORM’ 
nETURN 


igs 





PRINT 
ANSWER 


relish 
ieee 
OVL 
COM 
10 
SED Gat 
AAA 
BAA 
CAA 
AAT 
BAD 
ATT 
BDD 
ATA 
Jiakok 
NEE 
REL 
REL 
REL 
REL 
COM 


Programs Required for FL/1-&% 


> ADI 


APPENITIX B 


PL/1-8@ PROGRAMS 


DCLFUN 
UN 


TA 
PRINT 
ANSWER 


TA 


MONITOR 


Jeb ure 


PCLSUN 
TCLFUN 


» DCLFUN 


[TCLFUN 
CUCLFUN 
ICLSUN 
LCLFUN 


LATAINFO 


PLI 
Ee 
OVL 
REL 
eh 
AAA 
EAA 
CAA 
AAT 
BAL 
eeHD 
fe JB AE 
ADA 
TAT 
JEtal 
ESE 
KEL 
hob 
le. Bk 
SM 
Pi 


AA 
MONITOR 


TCLINP 


Le ph 
INPUT 
INPOT 
INPUT 
INPUT 
NGA) 
INPUT 


Operation 





MONITOR: proc options(main,stack(512) ); 


acl 
ftc entry (float binary) returns(char(17) var); 
ace i 
TIME fixed ext static init(2) 
M@lmOCK fixed ext static init( 
ACLOCK fixed ext static init ( 
ALT(8) float ext static init(o, 
DA(&) float ext static init(@ 
dummy _ptr ptr, 
Goce fixed(7), 
St lodt ext, /* correction factor for error signal 
datafile file; 


! 


aie 1 
Sinclude “dclmon.txt’; 


’ 


open file(datafile) record sequential output title(’data.flt’) 
env(b(512)); 
Gm = ftclexp(char(-delay/ALIMIT))); 
do while (DCLOCK <= DLIMIT); 
eo mireaiok LELOCK, LCLOCK’); /* ®BLIT */ 
code = 1; 
ACLOCK = @; 
panclude “ad.ttt’; 
Paraclude “dd.ttt’; 
code = 1; 
fomwolle (ACLOCK <= ALIMIT —- 1); 
meee Mee lot (ACLOGK, AGLOCE’ ,TIME, TIME’); /* ELI */ 
if ACLOCK = delay then do; 
winiciude “da.ttt ; 
end; /~ if */ 
Pyaclude “aa.ttt “; 
write file(datafile) from(dat); 
fein = LIME + 1; 
MopoGk = ACLOCK + 1; 
ena; /* do */ 
Meeock = ECGLOCK + 1; 
Smee /~ do */ 
e@leose file({datafile); 
end MONITOR; 


LOD 





PAN PL] */ 
mimelnude dclfun.ttt’; 
de |. 

ftc entry (float binary) returns(char(17) var); 
dic 1 

ioMcatiged..ext, 

Mm float ext, 

miatlioat ext, 

mio. o) variable entry(fixed(7)) returns(ptr), 

EA fixed static init(d), 

Mmotetixed Static init(¢), 

DA(8) float ext, 

ACLOCK fixed ext, 

code fixed(7), 

result ptr, 

m(5) label, 

MaweD oreo Dita C ptr,id ptr,ie.ptr) ptr; 
Zinclude “dclinp.txrt’; 


memelude “input.ttt’; 
mercode = 1 then do; 
mira < TIME then do; 
i meonririrs) Oo) & (nr inputs< €) 
then do; 
m(1)=la; 
m2 =p; 
m(3)=lc; 
m(4)=ld; 
m(5)=le; 
goto m(nr_inputs); 
end; /* then */ 
else goto lf; 


Mees ie ptr = input (5) (cod(5)); 
ld: id_ ptr = input (4) (cod(4)); 
eric o ptr — input (35) (cod(3)); 
Moeeioeptur = input (2) (cod {z)); 
taeda ptr = input (1) (cod(i)); 
lf: O =F; 

BAL = TIME; 


mimoelwde func.txt : 
’ 


KA = TIME; 


emda; /* if BA */ 
Smose /* if code then */ 
else 

if pbAl < TIME then do; 


O= F; 

BAL = TIME; 
ema, /* HAL */ 
Bereiude “out.ttt’; 
result = addr(F); 
return (result); 

end; 





Ve RY .PLI */ 
Minclude “dclfun.ttt’; 


are 1 

ftc entry (float binary) returns(char(17) var); 
dic | 

MA fired static init(d), 

T float ext, 

mi float ext, 

input(5) variable entry(fixed(7)) returns(ptr), 

CF float ext, 

mS) float eet, 

sode fj rea), 

Result ptr, 

m(S) label, 

Seep voro ptr,ic ptr,id_ptr,ie ptr) ptr: 


Mmaciude “dclinp.txt’; 
fic lude “input.ttt’; 


if EA > @ then do; 
mernpeinples> 0) & (nr_inputs< 6) 

then do; 
m(1)=la; 
mice =I; 
m(3)=1c; 
m(4)=ld; 
m(5)=le; 
goto m(nr_inputs); 

atid 7 then. */ 

else goto lf; 
mee eiemogr — input (5) (cod(5)); 
tmrcmpir =simput (4) (cod (4) ); 
Newmucmpitr = 1aput (3) (cod(s) ); 
Pi woe phr = input (2) (cod (2) ); 
Pimiaoor —minput (1) (eed (1 )); 
th ae 

pemeim@e fyunc.txt +; /* AL(X) = [A 5 BTC *, 


emg, /* then */ 
Pesuilt = addr(F); 
EA = 1; 


pemcitide out.ttt’; 
3 

return (result); 
end; 


Lore 





Poon PLL =/ 
Punelude dclfun.ttt’; 
aque 1 

pve ventvmy (float binary) returns(char(17) var); 
dic 1 

TIME fixed ext, 

T float ext, 

ieeetloat ext, 

input(5) variable entry(fixed(7)) returns(ptr), 

ED fixed static init(@), 

EDL fixed static init(@), 

CF float ext, 

AD(8) float ext, 

DCLOCK fixed ext, 

code fixed(7), 

Result ptr, 

m(5) label, 

eeu oe ptir,ic ptr,id ptr,ie ptr) ptr; 
%include “dclinp.txt’; 


meeclude “input.ttt’; 
if code = 1 then do; 
meet < ECLOCK then doa; 
Poor immpwmts> 9) & (nr_inputs< 6) 
then do; 
m(1)=la; 
ce) — 1), 
m(3)=lc; 
m(4)=ld; 
m(5)=le; 
Zoro minr inputs ); 
end; /* then */ 
else goto lf; 


Mere vemntr = input (5) (cod({5)); 
idee Td optr = input (4) (cod(4)); 
Neeeteeptc = input (3) (cod(3)); 
ieeeebeeor = input (2) (cod (2)); 
Maca laepth = input (1) (cod(1)); 
Bes O = fF; 


Bet = LOCLOCK; 
minclu@e “func.txt’; 


9 
tie = DCLOCK; 
ae, /* if EL */ 
end; /* if code then */ 
else do; 
mee eri < LCLOCK tnen do; 
= F; 
met = LOLOCS; 
amd; /* EL1 */ 
end; /* else */ 
me lude “out.ttt’; 
Gesult = addr(?F); 
HMecurn (result); 
end ; 


eS 





VOR PLT */ 
include “dclfun.ttt’; 


die | 

ftc entry (float binary) returns(char(17) var); 
del 

TIME fixed ext, 

tT float ext, 

Plo float ext, 

input(5) variable entry(fixed(7)) returns(ptr), 

Gir loat ext, 

mes) float ext, 

TA(8) float ext, 

code fixed(7), 

mesult ptr, 

m(5) label, 

Sreapeawoeptr,ic ptr,id ptr,ie ptr) ptr; 


memolude dclinp.txt’; 
Iaeoelude “input.ttt’; 


meet iMes > 8 then do; 
moeminr inputs> @) & (nr_inputs< 6) 

then do; 

m(i)=la; 

meece =) Ds 

m(3)=lc; 

m(4)=ld; 

m(5)=le; 

goto m(nr_inputs); 
end; /* then */ 
else goto lf; 


Wek iemptr = input (Sy) (1); 
ir idsptr = input (4) (1); 
Meeecm ptr — input (3S) (1); 
hosel omeie = input (2) (1); 
Mase len pur = input {1) (1); 
10 


Acme nimemext ; /* PACX) = {8B = ECT 


b 
ead, /* then */ 
meet = addr(r); 


memclude “out.ttt~; 


b] 
Return (result); 
end; 


MONS, 





REFORM: proc options(main); 


che d 
(onputeoutput) file, 
puff char(1) varying; 


open file(input) stream env(b(1024)) title(%$1.$17); 
open file(output) stream output env(bd(1024)) title("$2. 


do while (°1°b)3 
read file (input) into _(outt); 
purt =mieams late (buff, Nos 
write file (output) on (buff); 
end; 
end reform; 


PRINT: proc options(main); 


gic | 
mana dec(i5,6), 
mek,m) fixed bin(7), 
maume,dat,temp str) char(128) var, 
einio2mput,output) file, 
peemciude datainfo.txt’; 
V] 


Mmeemeriie(input) record sequeitial title(°LATA.FLT”) 
env 0 (522) 

open file(output 
env(b(512) ); 


)$ 
) 


stream output title(“°$1.$17) 


do K = © to timel; 
mime = ~ 3 
temp str=“TIME="|;char(k); 


do m = i to length(temp_str); 


Pewemoste temp sscr,m,i) = “ © then 
time = timej;substr(temp str,m,1); 
end; — 
time = “ jm 4 


i; time; 
write file(output) from(time); 
Gow j= i to nr data; 
read file(input) into (data); 


g 
dat = “m j iiname(j)it°= ticharida 
Meter le(outout) from(dat); 
end; 
e@ay /* do k */ 
end; 


ra 
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Figure 15C 


Terminated Ramp Input With No Corrections 
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Response to Variable Step Inputs with end Without Corrections 
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Figure 1¢Ff 


kesponse to Variable Step Inputs with and Without Corrections 
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